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THERMAL EXPANSIVITY AND DENSITY ( OF INDIUM 
By Peter Hidnert and Mary Grace Blair 


ABSTRACT 
The linear thermal expansion of cast indium of high purity was investigated 


between —190° and +100°C and the cubical thermal expansion between 0° and 
50°C. The following equation was derived for the linear expansion of indium: 


L,= Lf{1+ (28.93¢ + 0.01342) 10-*], 


vhere Lo represents the length of the metal at 0°C, and L, the length at any 
temperature, t, within the range of the observations. The average coefficients 
of cubical expansion of indium were found to be 77X10-* and 101 10-* per 
legree centigrade for the ranges 0° to 25°C and 28° to 50°C, respectively. These 

tlicients of cubical expansion are not equal to three times the coefficients of 
inear expansion for the corresponding temperature ranges on account of the 
anisotropy of indium. 

The density of cast indium was found to be 7.281 g/cm* at 22.6°C. The 
density of indium calculated from the best available data on atomic weights and 
lattice constants is 7.40 g/em’. Observed densities at 0°, 25°, and 50°C are also 
given in this paper. 


CONTENTS 
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I. INTRODUCTION 


In 1869 Fizeau [1]! reported the value 41.710~* per degree centi- 
grade for the coefficient of linear thermal expansion of cast indium at 
40° C and the value 0.4238X10-* as the change per degree centi- 
grade in the coefficient, or rate, of expansion. He computed the value 
45.9X10-° per degree centigrade for the average coefficient of expan- 
sion of indium between 0° and 100° C. The work described in the 
present paper was undertaken in order to obtain data on the linear 
thermal expansion of high-purity indium over a wide range of tem- 
perature (—190° to +100° C). Data on the cubical thermal expan- 
sion and density of indium between 0° and 50° C were also obtained. 

[Industrial applications of indium have been indicated by Murray 

[2], Williams [3], Linford [4], Dyer [5], and Ludwick [6]. A compilation 
a some of the physical properties of indium was given by West- 

rook [7]. 


! Figures in brackets indicate the literature references at the end of this paper. 
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II. MATERIAL INVESTIGATED 


A rod of cast indium loaned to the National Bureau of Standards 
by the Indium Corporation of America, Utica, N. Y., was used for 
the determinations of thermal expansion and density. A spectro- 
chemical analysis ? was made on a piece cut from the rod of indium, 
Traces of Ag, Ca, Cr, Mg, and Tl, small amounts (less than 0.01 per- 
cent) of Al, Cd(?), Co(?), Cu, Mn, Si, and Sn, and larger amounts of 
Fe, Ni, and Pb, were reported. The results of a chemical analysis * 
were: Fe, 0.01 percent; Pb, 0.01 percent; and Ni, not detected. The 
indium content was estimated to be 99.9 percent. 

X-ray diffraction patterns * obtained on a piece cut from the rod 
of indium, after the second expansion test had been completed, indi- 
cated large crystals. A subsequent examination (see footnote 4) of 
an etched section of the piece of indium indicated that the average 
diameter of the crystals was 1 to 2 mm. 


III. APPARATUS 


Five linear expansion tests were made on the rod of indium (10 mm 
in diam). The length of the sample was 250 mm for tests 1 and 2, 
which were made from 20° to 100° C with the precision micrometric 
thermal-expansion apparatus described by Souder and Hidnert [s). 
The details of the specimen mounting for use in the furnace are illus- 
trated in figure 4 of their publication. 

Tests 3° to 5 were made at low temperatures with the fused- 


quartz-tube thermal-expansion apparatus © described by Hidnert and 
Sweeney [9]. Liquid air, solid carbon dioxide in a mixture of equal 
parts of carbon tetrachloride and chloroform, and ice and water 
were used for low temperature baths at about —190°, —78°, and 0° C, 
respectively. 


IV. RESULTS 
1. THERMAL EXPANSION 


Figures 1 and 2 show the linear expansion and contraction curves 
of the sample of indium on two cycles of heating and cooling between 
20° and 100° C and on three cycles of cooling and heating at low 
temperatures. Table 1 gives the coefficients of expansion and con- 
traction obtained from these curves. 


2 By Dorothy Grosch, of the National Bureau of Standards. 

* By R. K. Bell, of the National Bureau of Standards. 

‘By H. C. Vacher, of the National Bureau of Standards. 

5 Before test 3, the rod of indium was cut to a length of 200 mm. 

6 This apparatus was found to be unsatisfactory for determinations of the linear thermal expansion of 
indium above room temperature on account of flow of the sample caused by the combined weight of the foot 
of the dial indicator ond the movable fused-quartz rod of the apparatus. 
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Figure 1.—Linear thermal expansion of indium on heating from room temperature 
to 100° C and contraction on subsequeat cooling from 100° C to room tempera- 
ture. 


ach number adjacent to a point on the curve represents the number of observations at the temperature. 
The variations in these observations were too small to be indicated on this scale. 


‘ABLE 1.—Coefficients of linear expansion and contraction of 


a sample of indium 


Average coeificients of expansion and con- 
traction per degree centigrade 
Treatment ® 
—()° —80° to} 20° to 60° to 20° to 
to 0° C | +20°C; 60°C | 100°C | 100° C 


x10-* | «10-6 -6 | x10-6 | «10-8 | KIN 
,/ Heating. a | 29.0 30.3 29.7 
Cooling - has cast : : 30. 5 31.6 31 
bar \Heated to 100° C and cooled slowly. __. + an ¢ zs : ap 
, Cooling__| Reheated to 100° C and cooled slowly 


Cooling. Na ‘ SS ee a 
4) Heating s° ooled to —80° C and heated to +20° ¢ 


5 one \Cooled to —190° C and heated to +20° C 





Values from equation 2 25. 27. ; 30. 0 


* All heat treatments on the cast sample were incident to the tests on thermal expansion. 
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Figure 2.—Linear thermal contraction on cooling to low temperatures and expansion 
on subsequent heating to room temperature. 


Each number adjacent to a point on the curve represents the number of observations at the temperatul 
The variations in these observations were too small to be indicated on this scale. 
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From the average coefficients of linear expansion and contraction 
obtained in the five tests on indium (table 1), the following linear 
equation was derived by the method of least squares: 


a,= (28.93+-0.026,¢) 10-°, (1) 
where a; is the instantaneous coefficient of expansion (contraction) or 
the rate of expansion (contraction) at any temperature, ¢° C, between 
—135° and +80°C. The probable error of a,is +0.4X107*. The 
differences in the coefficients of expansion (table 1) and the magnitude 
of the probable error of a, indicate that this metal failed to give 
reproducible results on heating and cooling. 

The following second-degree equation was obtained by integration 
of equation 1: 
L,=L{1+ (28.93t-+0.013,t2) 10-4. (2) 


In this equation, Z) represents the length at 0° C and L, the length at 
any temperature, ¢° C, within the range of the observations. 

After the determinations on linear thermal expansion had been 
completed, the cubical expansion of the sample of indium was deter- 
mined between 0° and 50° C by weighing the sample in air at room 
temperature and in water at 0°, 25°, and 50°C. Table 2 gives the 
coeflicients of cubical expansion. These coefficients of cubical expan- 
sions ’ are not equal to three times the coefficients of linear expansion 
of the sample of indium on account of the anisotropy ° of indium. 


TABLE 2.—Coefficients of cubical expansion of indium 


| Average | 
| coefficient | 

Temperature | of cubical 
| 


expansion | 
per degree 
centigrade 


range 


oy xX 10-6 
0 to 25 77 
25 to 50 101 
0 to 50 &Y 


The coefficient of linear expansion of indium (30.0107*) at 40° C, 
computed from equation 1, is considerably less than the coefficient 
41.710-® reported by Fizeau [1]. The large divergence between 
these values may be ascribed to different orientations of the crystals 
in the two samples. 

Frevel and Ott [10] determined the coefficients of linear thermal 
expansion of crystal indium along the a-axis and the c-axis from 
X-ray measurements on the cell dimensions of 99.9-percent-pure 
indium. They reported 56X10-* and 13X10~° for the coefficients of 
linear expansion along the a-axis and c-axis, respectively, for the 
range ® from —17° to +9° C. From these coefficients of linear 
’ Except the coefficient of cubical expansion for the range 0° to 50°C. 

Crystal structure of indium is face-centered tetragonal [10, 11). 


’ This range of temperature was reported in a private communication by L. K. Frevel, The Dow Chemi- 
Co., Midland, Mich. 





432 Journal of Research of the National Bureau of Standards 


— 


expansion, the computed coefficient of cubical expansion is 125 10> 
for the range from —17° to +9° C. The coefficients of cubical ex- 
are less than the value calculated from the 
[10] at a lower range of temperature. 


pansion given in table 2 
results of Frevel and Ott 


2. DENSITY 


Before the determinations of thermal expansion were made, the 


density of the 250-mm rod of indium was found to be 7.281 g/cm’ at 
22.6° C. The density of indium was also calculated from the equation 


p—nm(tiNe, 

ac 
where D is the density, n is the number of atoms included within the 
unit cell of the metal, w is the atomic weight of the metal, // is the 
atomic weight of hydrogen, No is the Avogadro number, and a and ¢ 
are the lattice constants of the metal. If the following values " are 
substituted in this equation: 


n=4 

w=114.76 

H=1.0080 

No=6.023 X 10” 

a=4.58,< 1078 cm 

c=4.94, X 1078 cm, 
the computed density is 7.39, g/em*®. This value agrees within 2 per- 
cent of the observed density. Observed densities reported in the 
literature [14] range from about 7.12 to 7.42 g/cm’. 

After the completion of the thermal-expansion determinations, the 
density of the 200-mm rod of indium was calculated at several tem- 
peratures from the data obtained in the determinations of cubical 
expansion. The results are given in table 3. 


TABLE 3.—Density of indium 


Temperature Density 
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MEASUREMENTS OF ULTRAVIOLET SOLAR RADIATION 
IN WASHINGTON, 1936 TO 1942 


By W. W. Coblentz and R. Stair 


ABSTRACT 


Continuing earlier measurements of ultraviolet solar radiation (RP877), a 
summary is given of similar measurements made during the years 1936 to 1942, 
in Washington, D. C. Marked variations in the spectral quality and total in- 
tensity of ultraviolet solar radiation, attributable to variations in atmospheric 
turbidity, and, to a less extent, to well-known cyclonic and seasonal changes in 
ozone concentrations in the stratosphere, are recorded. On the clearest days 
the biologically effective ultraviolet of wavelengths shorter than 3132 A, incident 
directly from the sun and the surrounding sky at midday, ranges from about 75 
microwatts per square centimeter (uw/cm?) in midsummer to about 12 pw/em? 
in midwinter. 


CONTENTS 


II. Discussion of ultraviolet solar intensity data--_---- 
teferences and notes._..-...------ 


I. INTRODUCTION 


The evaluation of the biologically effective component of ultra- 
violet in sunlight is still in its infancy. The difficulties and uncer- 
tainties involved in present-day measurements of the ultraviolet com- 
ponent of sunlight are somewhat comparable with those experienced 
in the beginning of the evaluation of the total incident solar radiation 
of all wavelengths, half a century ago. 

For some years the present writers have been engaged in an investi- 
gation of instruments and methods of evaluating the biologically 
effective component of ultraviolet solar and sky radiation useful in 
heliotherapy. 

In a previous paper [1] ' measurements were given of the biologically 
effective ultraviolet solar radiation of wavelengths 3132 A and shorter, 
incident normally in Washington, during the years 1934 and 1935. 
The present paper summarizes similar measurements made on some 
of the clearest days during the years 1936 to 1942, inclusive. 

The hereindescribed measurements of ultraviolet solar intensities, 
incident normally upon the earth’s surface, were obtained principally 
for three purposes: 

(a) to secure data on the intensities, incident directly from the sun 
and from a small contiguous portion of the sky, that may be expected 
on relatively clear days at a sea-level station in midlatitude (which in- 


' Figures in brackets indicate the literature references and notes at the end of this paper. 
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formation is of important in biological and photochemical irradiation 
experiments using a small exposure aperture) ; 

(b) to secure data required in working out methods of evaluating 
the biologically effective component of ultraviolet solar radiation 
[1]—required particularly in one method of standardizing ultraviolet 
dosage-intensity meters [4, 11]; and 

(c) to determine whether, with the accuracy attainable in the radio- 
metric measurements, the reduction factors (F, G@) used in the calcu- 
lations of ultraviolet intensities, will require modifications with changes 
in the spectral quality of the biologically effective ultraviolet radiation, 
which varies with the ozone concentration in the stratosphere, which, 
in turn, varies with the latitude, the season, the barometric pressure 
(movements of masses of polar air) and, less systematically, with the 
sunspot cycle [5, 6, 18]. 


II. APPARATUS AND METHODS 


In previous publications, descriptions were given of the instruments 
used in making ultraviolet solar-intensity measurements and _ the 
methods of evaluating the biologically effective spectral component 
of wavelengths 3132 A and shorter [1, 2, 7]. It will therefore suffice to 
recall that the intensity measurements were made with a spectrally 
selective radiometer [2] consisting of Ti-photoelectric cells sensitive 
to wavelengths 3500 A and shorter, and suitable glass filters that 
exclude ultraviolet radiation of wavelengths 3150 A and shorter. 

The photoelectric cell in its mounting is shown in figure 1 of reference 
7]. The aperture exposed to the sun and surrounding sky is about 22°. 
As shown in previous papers, at high elevations the amount of ultra- 
violet sky radiation is small (p. 142 of ref. [7]). On the other hand, a‘ 
a sea-level station even a slightly hazy atmosphere greatly increases 
the amount of ultraviolet sky radiation; but the total ultraviolet 
radiation from the sun and from a small part of the surrounding sky, 
incident normal upon the photocell receiver, is greatly reduced—see 
the measurements of April 13 and 16, 1942, depicted in figure 10. 

Using a photoelectric cell and filter method [1, 2], the separation of 
the amount of biologically effective ultraviolet solar radiation of wave- 
lengths 3132 A and shorter from the total amount of all wavelengths 
measured by a particular photoelectric cell (having a threshold of 
response at 3200 to 3500 A), requires measurements of the total 
intensity of the radiation to which the photoelectric cell responds and 
the amount transmitted through the filters at various solar altitudes. 
A knowledge is required also of: (a) the spectral response of the photo- 
electric cell, (b) the ultraviolet spectral transmissions of the filters, (c) 
the ultraviolet spectral-energy distribution at different solar heights 
(air masses), and (d) calibration of the ultraviolet meter, in absolute 
value, against a standard of ultraviolet radiation [8]. 

At the beginning of these investigations [20], based upon the spectral 
range of the erythemic and antirachitic responses, the intense emis- 
sion line in the mercury arc at 3132 A was adopted as the long-wave- 
length limit of biological effectiveness. In the meantime, it has been 
established that the spectral range of the tanning reaction extends to 
longer wavelengths [9], upon which basis it seemed logical to place 
the limit of biological effectiveness at 3200 A [10]. Another reason 
for considering the adoption of a threshold of biological effectiveness 
at longer wavelengths was the development of photocells sensitive 
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only to wavelengths shorter than about 3200 A, for use in measuring 
solar and sky radiation on a horizontal plane [11]. If reproducible 
photoelectric cells with the same threshold of response were available, 
ultraviolet radiometry would be greatly simplified. Thus far, how- 
ever, 1t has been impossible to make photocells having the same 
threshold of response. Hence, in view of the well-established use of 
3132 A in specifying the long-wavelength limit of biological effective- 
ness, it seems desirable to continue the use of this wavelength limit. 
This, of course, requires the continuation of the use of complicated 
factors for reducing the measurements obtained with various photo- 
electric cells and filters. In the present paper (figs. 4 to 9, inclusive) 
the data are given for both spectral ranges—to 3132 A and to 3200 A, 
respectively. 

The procedure employed in calculating the factor P, used in cali- 
brating the photoelectric cell against a primary standard of ultra- 
violet radiation, and the factor G, used in evaluating the ultraviolet 
solar radiation of wavelengths shorter than 3132 A (or 3200 A) relative 
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FiguRE 1.—Reduction factors, G, used in determining the amount of ultraviolet solar 
radiation of wavelengths shorter than and including 3200 A relative to the total spectral 
range measured with photoelectric cells Ti-1 and Ti-4. 


to the total range of wavelengths to which the photocell responds, is 
described in previous papers [1, 13]. As a matter of record, in figure 1 
are illustrated the reduction factors G for the spectral band of 3200 A 
and shorter wavelengths for comparison with figure 5 of reference [1] 
for the band of 3132 A and shorter wavelengths. Obviously, if 
photocells can be developed having the same threshold response (say 
3200 A) for use as comparison standards, the same factor G can be 
used for all or be eliminated entirely by comparison with primary 
standards. 

The weak point in the calculation of the factor G is the lack of exact 
knowledge of the spectral distribution of energy in the extreme ultra- 
violet of solar radiation at the observing station. The method of 
deducing the solar spectral-energy distribution by means of various 
photoelectric cells and filters is described in previous publications 
((7] and p. 323 of ref. [1]). To simplify the calculations, a smooth 
spectral-energy distribution is assumed, although from the very be- 
ginning [20] it was noted that the energy curve is deeply indented 
(see fig. 2 of ref. [7] and p. 167 of RP1367 [18]) by numerous Fraun- 
hofer and ozone absorption lines. 





438 Journal of Research of the National Bureau of Standards 


This method of checking the ultraviolet spectral-energy distribution 
may be criticized because of the difficulty in determining the threshold 
and the shape of the spectral photoelectric response. The recent 
determinations of the spectral responses of our photoelectric cells, using 
a monochromator and a mercury- or cadmium-arc lamp, are therefore 
supplemented by integrated filter-transmission measurements of the 
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Figure 2.—Diurnal variation in transmission of ultraviolet solar radiation through 
standard filters, from January to June 1939, observed in Washington. 


See fig. 3 for seasonal variation. 


ultraviolet radiation emitted by a tungsten-ribbon-filament lamp 
operated at a high temperature (2,900°K), and using the spectral- 
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Ficure 3.—Diurnal and (dotted curves from fig. 2) seasonal variation in transmission 
of ultraviolet solar radiation through standard fillers, observed from July to Decem- 
ber 1989, in Washington, D. C. 


energy data observed by Hamaker [14] (communicated by Ornstein 
[14]) and by Hoffmann and Willenberg [15]. 

As shown elsewhere [16, 17], based upon integrated ultraviolet- 
transmission measurements made with four filters of known spectral 
ultraviolet transmissions and eight photoelectric cells differing widely 
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in spectral photoelectric response, the average values of the spectral 
emissivities of tungsten, in the band of wavelengths extending from 
2,400 to 3,500 A, observed by Hoffmann and Willenberg appear to be 
substantially correct. Hence, as mentioned above, in recent deter- 
minations of the spectral response of a photoelectric cell, the wide gaps 
between the response measurements, made with the spectrally isolated 
emission lines of the mercury and the cadmium are (at 3,024, 3,132, 
3,261, 3,342, and 3,663 A) are filled in on the basis of a smooth spectral 
photoelectric response curve that is adjusted to give calculated filter 
transmissions in agreement with the integrated filter transmissions ob- 
served on the ultraviolet radiation from a tungsten-ribbon lamp oper- 
ated at a known temperature. 

From our most recent studies it appears that the ultraviolet solar 
spectral-energy distributions, published in table 2 (p. 323 of ref. [1]), 
remain unchanged for air masses m=2, and less; but for air masses, 
m=2.3 and larger, the data appear to belong to air masses about 0.3 
larger, e. g., m=3.3 instead of m=3.0. However, for these large air 
masses the spectral intensities of the shortest wavelengths are so small 
that this uncertainty seems relatively negligible in the present-day 
stage of evaluation of biologically effective ultraviolet solar radiation. 
As noted in earlier papers [1, 12], the good agreement in the measure- 
ments of the biologically effective component of ultraviolet solar radia- 
tion of wavelengths shorter than 3,132 A, obtained with spectrally 
nonselective thermopiles and with wavelength-selective radiometers 
differing greatly in range of spectral photoelectric response, appears 
to be real and not the result of a fortuitous combination of factors, 
including (a) the spectral responses of the photoelectric cells, (b) the 
spectral transmissions of the filters, and (c) the transmissions of the 
integrated ultraviolet solar radiation through the various filters. 

In figures 4 to 10, inclusive, the ultraviolet intensities are given in 
the two above-mentioned spectral ranges, extending from the shortest 
incident wave lengths (at about 2900 A) to 3132 A and to 3200 A, 
respectively. The values of the intensities in these two spectral 
ranges differ by about 40 percent, varying from about 33 percent for 
air mass, m=8, to about 48 percent for m=1.05. In these illustra- 
tions the scale on the right-hand side is to the nearest rounded num- 
bers. Based upon some preliminary trials [4], it is believed that when 
satisfactory photoelectric cells are obtainable for use as primary 
standards, the calibration of meters for measuring ultraviolet solar 
and sky radiation incident on a horizontal plane will be greatly 
simplified. 

During the 7 years of measurements included in the present report, 
two of the filters (Ba-1 and Ba-3, a soft glass) became slightly 
scratched from frequent wiping with a fine-haired brush and were 
replaced by others of practically the same thickness. No difference 
was observed in the ultraviolet intensities deduced from the measure- 
ments with the old and new filters. Furthermore, frequent calibra- 
tion against a primary standard of ultraviolet radiation showed no 
systematic change in the sensitivity of the photocells (Ti-1 and Ti-4, 
illustrated in figs. 1 and 2 of ref. [7]) used throughout these researches. 
From this it appears that the instrumental precision from year to 
year was maintained within 10 percent, or even closer. 
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The accuracy attained in the measurements of ultraviolet solar 
intensities was limited entirely by meteorological conditions. As 
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Figure 4.—Measurements of ultraviolet solar radiation of wavelengths shorter than 
$182 A, also 3200 A and shorter, during some of the clearest days in 1936. 


noted in previous papers, on the clearest days at high elevations 
(Flagstaff, Ariz.) where the measurements could be repeated to 0.1 
to 0.2 percent, a change in spectral quality and total intensity, result- 
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Figure 5.—Measurements of ultraviolet solar radiation of wavelengths shorter than 
3132 A, also 3200 A and shorter, during some of the clearest days in 1987. 


ing from a change in solar height, was observed from minute to 
minute. 
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In contrast, at sea-level stations (San Juan, P. R., and Washington, 
Dp. C.) with a continual variation in turbidity in the lower atmos- 
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‘RE 6.— Measurements of ultraviolet solar radiation of wavelengths shorter than 
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phere, the precision attained is best shown graphically in the measure- 
ments of April 13 and 16, 1942, figure 10, where each measurement (the 
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instantaneous value in the group, say at m=1.2 and at m=1.3) is 
probably accurate to 1 percent. Nevertheless, because of atmos- 
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Ficure 8.— Measurements of ultraviolet solar radiation of wavelergths shorter than 
3182 A, also 3200 A and shorter, during some of the clearest days in 1940. 
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FiguRE 9.— Measurements of ultrqiolet solar radiation of wavelengths shorter than 
3182 A, also 3200 A and shorter, during some of the clearest days in 1941. 


pherie conditions the average values for these two days are best 
represented by the two lines drawn through these measurements. 
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Tbe wide spread in the values of the ultraviolet solar intensities 
depicted in figures 4 to 10, inclusive, is attributable to constantly 
varying meteorological conditions. This was shown in an earlier paper 
(fig. 6 of ref [1], where a detailed description is given of the effect of a 
slight change in turbidity upon the ultraviolet intensities observed 
throughout an exceptionally clear day, May 18, 1935. 

The accuracy attained in the absolute values of the herein-described 
measurements of ultraviolet solar intensities is somewhat of a con- 
jecture. As noted in a previous paper (p. 335 and fig. 6 of ref. [1]), 
using three distinctly different methods of evaluation (thermoelectric 
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Figure 10.— Measurements of ultraviolet solar radiation of wavelengths shorter than 
3132 A, dvring some of the clearest days in 1942. 


he large black dots, and the circles, connected with lines, indicate measurements on 2 days, differing 
greatly in atmospheric turbidity. 


and photoelectric), the ultraviolet solar intensities, of wavelengths 
3132 A and shorter, observed at midday in Washington during the 
clearest midsummer weather, range from 75 to 78 ww/cm?. 

Using measured intensities of ultraviolet solar radiation, from 
erythematogenic tests made at two sea-level stations (San Juan, 
P. R., and Washington, D. C.) and at two high elevations (7,300 and 
10,000 ft.; Flagstaff and the San Francisco Peaks, Ariz.), the uncer- 
tainty in the radiometric evaluation seems to be no greater than the 
uncertainty in the biological evaluation used in checking the radio- 
metric measurements, viz., 10 to 15 percent [1]. This uncertainty 
is due to a lack of knowledge of the spectral-energy distribution of 
solar radiation in the biologically effective range of wavelengths 
3200 A and shorter. Similar radiometric and biologic comparisons 
made by Greider and Downs, by Luckiesh, and by Pettit (discussed 
on p. 973 of ref. [20], reduced to the same range of wavelengths, are 
in good agreement with the foregoing data. 


III. DISCUSSION OF ULTRAVIOLET SOLAR INTENSITY 
DATA 


In previous papers it was shown that, for the same solar altitudes, 
the atmosphere is more transparent to ultraviolet radiation of short 
wavelengths in the afternoon than in the forenoon and more trans- 
parent in the autumn than in the spring (p. 144 and fig. 5 of ref. [7]). 
The reason therefor is the presence of less ozone in the stratosphere in 
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the autumn than in the spring [5, 6], and presumably less in the after- 
noon than in the forenoon. While, in the present stage of evaluation of 
ultraviolet solar radiation, this difference in atmospheric transparency 
is not taken into consideration (see figs. 2 and 5 of ref. {1]), it is evident 
that, at stations where atmospheric conditions are sufficiently con- 
stant to warrant such refinements in reducing the observations, this 
variation in ultraviolet transparency should be included in the 
calculations. 

In figure 2 of a previous publication [1] a sudden, presumably sea- 
sonal, change in atmospheric transparency was observed in the tropics, 
Similar changes are evident from the filter measurements made at 
Flagstaff, Ariz., elevation about 7,300 ft (fig. 5 of ref. [7] and fig. 3 of 
ref. [18]). In spite of the great variation in the dirunal measurements, 
the average values of the filter transmissions observed from January 
to June 1939 (fig. 2 and the dotted curves in fig. 3) are conspicuously 
higher than those observed from July to December 1939 (fig. 3). 
This seasonal difference in filter transmissions is equivalent to a differ- 
ence of about 0.2 air mass (3° to 5° in solar height) or a decrease of 
perhaps 8 to 10 percent in the amount of ozone in the stratosphere in 
the autumn. 

The measurements of ultraviolet solar radiation of wavelengths 
shorter than and including 3132 A, also in the spectral range of 3200 A 
and shorter wavelengths, made on some of the clearest days during the 
years 1936 to 1942, inclusive, in Washington, are depicted in figures 4 
to 10. 

The wide difference in the ultraviolet intensities observed from day 
to day at any given solar height (air mass in these illustrations) is 
attributed primarily to variations in atmospheric turbidity, although 
changes in the amount of ozone in the stratosphre also cause varia- 
tions in ultraviolet intensities. Even a slight increase in atmospheric 
scattering (as judged by the solar corona) greatly decreases the 
intensity of ultraviolet radiation incident on a surface exposed 
normal to the sun’s rays. This was repeatedly observed in the meas- 
urements of 1942 (fig. 10), the most of which were made in connection 
with the calibration of a photoelectric cell used for measuring ultra- 
violet solar and sky radiation incident on a horizontal plane (4, 11]. 
In the course of these measurements of solar and sky radiation it was 
repeatedly observed (by shadowing the receiver from the direct solar 
rays) that it was only during the noon hours of the clearest days in 
early summer that the ultraviolet radiation from the whole sky inci- 
dent on a horizontal surface was slightly less than the directly 
incident solar rays, and that when the sky was slightly turbid, the 
ultraviolet radiation from the whole sky was at least 2.5 times greater 
than the directly incident solar rays [3, 4]. 

A specific example of the decrease in intensity caused by atmospheric 
scattering is illustrated in the measurements of April 13 and 16, 1942, 
depicted in figure 10. On April 13 (sky blue; no corona) the ratio of 
the ultraviolet in sky radiation to direct solar radiation was 1:1. On 
April 16 (sky turbid; white corona) the ultraviolet sky radiation was 
2.5 times greater than the direct solar radiation. In this instance, 
however, contrary to general experience, the filter transmissions were 
lower on April 16 than on April 13, indicating the presence of a rela- 
tively greater amount of ultraviolet of the shortest wavelengths on 
April 16. This might be due to selective atmospheric (Rayleigh) 
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scattering (see p. 146, of ref. [7]). On the other hand, the barometer 
was higher (760 mm for 2 days; hence, there probably was less ozone) 
on April 16 than on April 13, when the barometer (759 mm) was rising 
from a (cyclonic) barometric low that occurred on April 10 to 11. It 
is therefore probable that on April 16 the reduction in ultraviolet 
intensities by atmospheric scattering was compensated for appreciably 
by an increased transparency (caused by a decreased ozone content) 
of the stratosphere to ultraviolet solar radiation. In this connection, 
reference is made to measurements procured at high altitudes (Flag- 
staff, Ariz., elevation 7,300 ft) showing the effect of (a) scattering by 
atmospheric water vapor and (b) variation in barometric pressure, 
and an accompanying variation in ozone content of the stratosphere, 
upon the incident ultraviolet solar radiation (see figs. 1 to 5, inclusive, 
of ref. [18]). 

With constantly changing atmospheric pollution and cloudiness, 
especially at sea-level stations, great variations in ultraviolet inten- 
sities may be expected. Hence, the measurements depicted in figures 
4 to 10 represent only what may be expected under the most favorable 
atmospheric conditions. 

As already mentioned, because of the seasonal decrease in ozone 
content and, hence, increase in transparency of the stratosphere, the 
ultraviolet intensities should be systematically higher in the autumn 
than for similar solar heights at other seasons. This is shown in 
figure 4 (Sept. 25, Oct. 20), figure 5 (Nov. 9, Dec. 3), figure 7 (Sept. 
5, Oct. 16), figure 9 (Oct. 8, Nov. 17), ete. However, while high 
ultraviolet solar intensities seem to predominate in the autumn, they 
occur frequently in the late winter and early spring (e. g., February 
and April; figs. 5, 6, 7, and 8), when the highest concentration of ozone 
usually prevails in the stratosphere (see fig. 1 of ref. [4]). Correcting 
for solar distance, the ultraviolet intensities in March, relative to 
September, would be reduced by about 2 percent. Evidently this 
does not explain the high ultraviolet solar intensities in February, 
which, more logically, are explainable on the basis of a temporary 
reduction in the ozone content (the ‘‘ozone holes” observed by 
Kiepenbeuer [19]) and a concomitant increase in transparency of the 
stratosphere to ultraviolet radiation of very short wavelengths [19]. 

These measurements were begun when the relative numbers of 
sunspots were low (minimum in 1934). They include the maximum 
of the sunspot cycle (April 1937; secondary maximum January 1939) 
and the period of rapidly decreasing sunspot numbers in 1941 and 
1942—predicated minimum about 1944. 

The clearest weather is usually associated with a high or average 
barometric pressure. It is therefore rare to secure ultraviolet solar- 
intensity measurements during a cyclonic disturbance, when atmos- 
pheric ozone is higher than normal (see fig. 5 and discussion on p. 575 
to 577 of ref. [18]. 

The seasonal variation in atmospheric ozone, especially in midlati- 
tude (see fig. 1 of ref. [4]), also must be considered. Since the time 
of making the hereindescribed measurements, from season to season, 
could not be selected (for example, there were but few clear days 
during the summer of 1940—see fig. 8, m=1.05) only generalized com- 
parisons can be made. It therefore seems illusory to attempt to 
seek a correlation of these data with the sunspot cycle. For example, 
in the measurements of 1938 (fig. 6) the noon-hour ultraviolet inten- 
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sities (u.v.Q.) for m=1.1 are higher than form=1.05. This was owing 
to clearer weather in April than in June, and not because of the effect 
of sunspots upon the ultraviolet intensity. Selecting 2 years having 
the same sunspot numbers, in 1936 the u.v.Q. for m=1.05 is higher than 
in 1941; but in September the intensities are the same as, or they are 
higher than, observed during some other years. Again, in comparison 
with 1941, the u.v.Q. values of 1937 are low in April and September: 
and they are practically the same in February, November, and 
December. 

In the absence of factors for eliminating the effect of atmospheric 
turbidity and the variability (seasonal and meteorological) of atmos- 
pheric ozone, it is impracticable to determine from these data the 
elect of sunspots upon the amount of ultraviolet emitted. In the 
few instances when observations could be made near the time when 
the ionosphere (and radio transmission) was disturbed by large out- 
bursts of sunspots, the ultraviolet intensities showed nothing con- 
spicuously different from the usual values. 

The measurements depicted in figures 4 to 10 indicate what inten- 
sities may be expected in irradiation experiments using biological 
effective ultraviolet radiation incident from the sun and from a narrow 
portion of the surrounding sky, practically free from pollution by 
smoke, on relatively cloud-free days, at a practically sea-level station, 
Washington, D. C. The great variations in intensities from day to 
day, for the same solar height (air mass, in these illustrations), is 
caused primarily by differences in atmospheric turbidity. <A specific 
example is illustrated in figure 10; April 13 and 16, 1942, described on 
a preceding page. At the same solar height (air mass traversed) the 
intensity differs by 30 to 35 percent, primarily because of atmospheric 
scattering. As will be shown in a subsequent paper, the total incident 
ultraviolet solar and sky radiation, integrated for the whole day, 
was 16 percent less on April 16 when the sky was more turbid. 

In figures 4 to 10, inclusive, the dates attached to the low intensities 
indicate days that were conspicuously hazy; with no attempt at cor- 
relating the low values with atmospheric scattering and with in- 
creased ozone in the atmosphere as indicated by a decrease in baro- 
metric pressure. 

From the data depicted in figures 4 to 10, inclusive, it appears that 
at a sea-level station in midlatitude, on exceptionally clear days, the 
noon-hour intensity of biologically effective ultraviolet radiation of 
wavelengths less than 3132 A, incident directly from the sun and a 
small portion of the surrounding sky, varies from about 75 yw/cm’ 
in June (air mass, m=1.05) to about 12 ww/cm? in December (m= 
2.1). For similar measurements, made during 1934 and 1935, with 
noon-hour designations of intensities for all the months, reference is 
made to figure 7 of ref. [1]. 

In a previous paper (fig. 6 of ref. [4]) data were given showing that, 
in high latitudes (78° N) at the noon-hour on the clearest days in mid- 
summer, the intensity of the biologically effective ultraviolet radiation 
of wavelengths shorter than 3200 A from the sun and the whole sky, 
incident on a horizontal plane, ranged from 30 to 40 ww/cm*, As 
shown in the figures 4 to 10.of the present paper for the band of wave- 
lengths shorter than 3132 A, the corresponding values range from 10 
to 14 uw/em?. The longer spectral range includes more of the ultra- 
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violet rays that produce tanning. The shorter spectral range is of 
more significance in the prevention of rickets. 

As noted in the previous paper [4], the intensity and biological 
effectiveness of ultraviolet solar and sky radiation, in midsummer, in 
high latitudes (78° N) is not markedly different from that observed 
at the same solar height in December (m=2.1) in Washington. It 
was shown also that while these intensities are low, if long exposures 
can be made, they are effective in heliotherapy. 
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A TUNGSTEN.IN-QUARTZ LAMP AND ITS APPLICATIONS 
IN PHOTOELECTRIC RADIOMETRY 


By R. Stair and W. O. Smith* 


ABSTRACT 


This paper gives constructional data relating to a tungsten-filament-in-quartz 
lamp and a discussion of its applications to precision filter radiometry. This 
lamp is adapted for use in the calibration of the spectral and total radiation sensi- 
tivity of phototubes. In its applications, it supplements the spectroradiometer 
for calibrating phototubes while employing the tubes with the associated filters 
used in measuring the radiation under investigation. 

The spectral radiation output from such a lamp depends upon the temperature 
of the filament, the emissivity of tungsten, and the transmission of the fused- 
quartz envelope. A table of relative blackbody radiation intensities has been 
calculated for the temperature range 2,500° to 2,900° K in the wavelength 
interval 2300 to 3500 angstroms. 
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I. INTRODUCTION 


The phototube has attained a place of great importance for making 
radiation measurements in the ultraviolet, visible, and near-infrared 
regions of the spectrum. The field of application is continually 
extending. Electronic devices for the amplification of the photo- 
electric current allow measurements of energies many thousand times 
smaller than are possible with the most sensitive thermopiles and gal- 
vanometers available. Reliable results are obtained by using stable 
amplifiers with microammeters or stable d’Arsonval galvanometers in 
cases where, because of mechanical vibrations, temperature changes, 
air currents, and the like, it is impossible to employ a sensitive thermo- 
pile-galvanometer radiometer. 

Phototubes, however, have a selective response with respect to 
wavelength, which varies from tube to tube and also from point to 
point on the surface of the same electrode, both as to absolute inten- 
sity and as to wavelength range of sensitivity, even in tubes of the 
same type and process of manufacture. It is, therefore, necessary not 
only to calibrate individual tubes for the specific purpose to which they 
are to be applied, but the identical portion of the surface of the photo- 

ibe that is to be used must be calibrated. 


* Associate Physicist, Bureau of Plant Industry, U. 8. Department of Agriculture. 
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None of the existing procedures for the calibration of phototubes 
appears to be adequate. Accurate calibration has become an impor- 
tant problem in the field of photoelectric radiometry. Radiation from 
an ordinary Mazda CX lamp, used with filters, has been employed in 
the spectral calibration of some of the phototubes used for the measure- 
ment of short wavelength ultraviolet solar radiation and the deter- 
mination of the amount of ozone in the stratosphere. [1].! More 
recently a tungsten-ribbon-filament lamp, with a fused-quartz window, 
has been used in the measurement of the spectral transmission of 
spectroradiometers [2]. The limitations of transmission of glass are 
removed, but since auxiliary diaphragms must be employed, in many 
cases it is possible to irradiate only a smal] part of the active surface of 
the phototube. The sensitivity of different parts of the active surface 
is known to vary considerably, and calibrations made for a portion of 
the active surface are not valid when the entire, or even another, por- 
tion of the active surface is used. 

To eliminate the difficulties just mentioned, a tungsten-filament 
lamp, enclosed by a fused-quartz envelope, was constructed. Since 
this type of lamp appears to be a satisfactory standard source of con- 
tinuous radiation, it will be fully described, and its properties and use 
discussed. Other types of tungsten-filament lamps with fused-quartz 
envelopes have been designed and used for other purposes [5, 17] with 
practical results. 


II. DESCRIPTION OF LAMP 


The lamp constructed is shown in figures 1 and 2.? It consists 
of an 8-mil straight filament of pure tungsten wire enclosed in a spher- 
ical fused-quartz envelope 4 inches in diameter. A straight round 
wire rather than a coiled filament was chosen to insure relatively 
even temperatures at all points on the wire. Black tungsten shields 
cover the filament ends and loops, as well as the molybdenum sup- 
porting hooks, and prevent radiation from these cooler areas of the 
filament reaching the receiver being calibrated. The filament con- 
sists of four evenly spaced “hairpin” loops arranged, in a plane, over 
an area approximately % by 1% inches. The blackened tungsten 
shields were placed inside the lamp for simplicity in the use of the 
lamp. It was realized that these shields would become heated and 
reradiate energy at long wavelengths. Since not more than 10 
percent of the total energy radiated from the lamp filament can be 
absorbed by each shield, and since the surface of each shield is more 
than 10 times the total surface area of the filament, it follows that the 
total energy radiated per unit area of the shields is less than 1 per- 
cent of that for the filament. Hence the temperature reached by the 
shields cannot be very high, and therefore the spectral distribution 
of the radiation is limited to long wavelengths in the infrared, which 
do not excite available phototubes. 

The lamp is filled with nitrogen at a pressure of 0.6 atmosphere. 
The filament was “flashed” in nitrogen, prior to sealing it in the 
fused-quartz envelope, at a temperature of about 3,000° K for 1 hour 

! Figures in brackets indicate the literature references at the end of this paper. 


2 The advice and assistance of R. A. Nelson, of the Bureau of Plant Industry, in the construction of this 
lamp is gratefully acknowledged. 
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Figure 1.—Tungsten-in-quartz lamp, 
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to free“it of “any impurities. It has been shown [2, 5] that different 
lots of tungsten agree in spectral emissivity, probably becausefat 
Ligh temperatures most impurities are evaporated, leaving practically 
pure tungsten [5, 7]. To extend the useful life of such a lamp, it 
should not be operated at a temperature exceeding 2,800° K after 
flashing. 

A lamp of this type, as shown by experience, should have a power 
input of about 500 watts to insure sufficient radiation output for the 
calibration of the relative spectral response of phototubes by the filter 
method, especially the types of tubes sensitive only to ultraviolet 
radiation. The problem of a source of power and its control must be 
considered when designing a lamp. It is convenient to have a lamp 
operating somewhat below 100 volts and at the same time requiring 
less than 10 amperes, so that it may be controlled with the common 
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FiguRE 2.—Constructional details of the lamp. 
A, Shields; F, filament; Q, quartz supports. 


type of laboratory equipment. This is realized by using about 16 to 
18 inches of 8-mil tungsten wire as a filament. A short ribbon filament 
operated at a low voltage and high current requires special control 
equipment. 

The entire lamp envelope, except for the Pyrex press, is constructed 
of high quality fused quartz, since it has a high transmission through- 
out the region from below 2500 to above 20000 A. This construc- 
tion, together with shields which prevent radiation from the cooler 
ends and loops of the lamp filament reaching the phototube, eliminates 
the need of an additional diaphragm, or lens, between the lamp and 
tube under examination and thus the phototube is completely and 
uniformly irradiated with radiation of known spectral quality. In- 
homogeneity and variations in the thickness of the quartz envelope 
render an image of the lamp filament distorted (see Fig. 1). This 
makes no difference, however, in the use of the lamp in the calibration 
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of phototubes. If it is desired to use an image of the filament for 
some purpose, the lamp could be provided with a sealed-in window 
of polished fused quartz. 

Finally, this lamp may be operated at any temperature from 2,000° 
K to about 2,800° K without excessive depreciation of the tungsten 
filament or blackening of the lamp envelope. 


III. SPECTRAL-ENERGY DISTRIBUTION OF THE LAMP 


In applications such as the calibration of phototubes, an accurate 
knowledge of the distribution of the energy radiated from the lamp 
is require ed. For the particular type of lamp considered in this paper, 
it is a function of the filament temperature, the emissivity of the mate- 
rial composing the filament, and the transmission of the envelope. 
For this purpose, therefore, blackbody radiation, the emissivity of 
tungsten, and the transmission of the fused-quartz envelope must be 
considered. 

The spectral quality of the radiation from a blackbody depends 
only upon the temperature. The intensity J, of the energy at wave- 
length d, according to Planck’s radiation law, is given by 


Jy=O (eC T—1)7, 


where C, is the first radiation constant, C, is the second radiation 
constant, 7'is the absolute temperature, and ¢ is the base of the natural 
logarithms. If \ is in microns, and 7’ in degrees Kelvin, then C, has 
the value 3.732 10* watt cm? and C, the value 14,3604 deg [3]; J 
is then the radiant power, in watts, over a solid angle 27, for unit area 
(em?) of the blackbody per micron of wavelength. Values of J, 
can be calculated directly from eq 1. The value of the radiation 
constant C, was taken as 1.436 em deg (14,360 uw deg), which has been 
shown in a recent summary by Wensel [3] to be the most probable 
value. 

Values of the relative intensities, from eq 1, for wavelengths 2300 
to 3500 A and for temperatures of 2,500° to 2,900° K are given in 
table 1. This table was calculated for use in studies of the short 
wavelengths of ultraviolet solar radiation and should be very 
helpful when ultraviolet radiation from the lamp is used.’ Values 
of J, relative to that for \=0.3500 micron are given for the tempera- 
ture range 2,500° to 2,900° K in 50° increments. Wavelengths are 
given from 2300 to 3500 A in increments of 50 A. 

The values given in table 1 are (J,/Jo.3500) X (10%); or what is 
equivalent, Jo.3500 has been set equal to 10,000. The constant C,, 
it is to be observed, is not used in the calculations, since these are rela- 
tive values. Values of ¢* used to compute the data of these tables 
were taken from tables [11] recently published by the Work Projects 
Administration. 

The values tabulated in table 1 are, in most cases, considered to 
be accurate to better than one in the final ee figure. 


Complete tables to 4 significant figures over the range 2,000° to 3,100° K, in 25° steps, in the wavelength 
interval 2300 to 3500 angstroms, may be obtained upon application to the Nation ul Bureau of Standards. 
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[aBLE 1.—Relative values of J, based on a value of J, at X=0.3500 micron of 10,000 


Blackbody temperatures 


Wave 
lengths -- al 5 . 
() | 2,500° K | 2,550° K | 2,600° K | 2,650° K | 2,700° K | 2,750° K | 2,800° K | 2,850° K | 2,900° K 


22 25 
32 38 
48 55 
ao 
86 as 110 


122 136 152 169 | 188 | y | 228 


169 | 188 208 230 253 | y 303 
231 =| 255 280 307 336 | 366 397 


| 311 341 373 | 406 440 | 5 | 514 

380 | 415 | 451 490 | 529 571 | 614 | 659 

504 547 591 636 | 683 732 =| 3 | 835 
662 | 712 | 764 818 873 | 930 9} 1048 | 1109 
859 918 | 978 1040 1104 | 1169 235 | 1303 1371 
1103 | 1171 | 1240 1311 | 1383 1456 5x 1605 1681 
1402 | 1480 1558 | 1637 1716 1798 1962 2045 
3050 | 1767 1853 | 1940 | 2027 | 2114 2202 22h 2379 =| 2468 
3100 | 2207 2301 2395 | 2488 2582 2676 | 2770 | 2863 | 2956 
3150 2734 | 2834 295: 3032 3130 3227 3324 3421 3516 
3200 3361 3464 | 3565 3666 3766 3865 3963 | 4059 | 4155 
3250 4099 4202 303 4403 4501 4597 9: 4786 4878 
3300 4964 5062 | 5157 4250 5343 5434 552% 5609 5604 
3350 5970 6057 | 6141 6224 6304 6382 | 45¢ 6534 6607 
3400 7134 7201 7267 7332 7393 7454 51% 7569 7625 
3450 $472 | 8511 8549 8586 8622 8657 8723 8754 
3500 10000 10000 10000 10000 10000 10000 10000 10000 10000 


Data [12, 13, 14, 15] are to be found elsewhere that contain similar 
blackbody calculations for longer wavelengths. The tables of Skog- 
land [13], for example, list values of (J)/Jo.s900 from 3200 to 7600 A 
for a temperature range of 2,000° K to 3,120° K. A slightly different 
value of the second radiation constant C, is used in most published 
work, but a correction factor has been determined [14]. Where it is 
desired to calculate a specific value, a short method of calculation given 
by Coblentz [12] is useful. 

The emissivity of tungsten has been determined by several workers 
(4, 5, 6, 7] for different temperatures. The data which appear to best 
represent the published values for the temperature range 2,700° to 
2,900° K are tabulated in table 2. These values may be used in 
work of this kind throughout the range 2,500° to 2,900° K, since 
only relative values enter into the calculations. 

As the transmission in the spectral region 2500 to 20000 A of 
fused quartz of the thickness (about 1 mm) employed in this lamp 
is nearly uniform, no correction is necessary for most practical pur- 
poses. The transmission of the fused-quartz envelope of the lamp 
must be taken into consideration, however, when a phototube has 
its major response in the short wavelengths, for example, a tungsten 
or tantalum phototube in a quartz or special ultraviolet transmitting 


glass envelope. 
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TABLE 2.—-An example of the application of the tungsten-in-quartz lamp in th: 
determination_of the_spectral response of a phototube 


(The spectral-response data (col. 5) are adjusted until the observed and calculated filter transmissions 
are in close agreement] 








Rela- | Relative radiation fro y 
tive tungsten times relatiy; 
radia- response of phototube 
te tion times filter transmissio; 
Rela- Boga ; - _| from 
tive | re iv | tung- 
radia- | Emis- , sten 
tion | sivity : s | times 
from | factor | ~ f rela- 
black | for ee tive 
| body tung: | 5 so . | | re- 
| > yhoto-} np, : in, a.2 | Sponse 
be ‘ ie | Ba-l | Ba-3 | of 
A | | photo- 
tube 
(col. 
4Xcol. 
5) 


Filter transmissions 














l 


A | 
2560 | 
2580 | 
2600 | 
2620 
2640 | 
2660 
2680 
2700 
2720 
2740 
2760 
2780 
2800 
2820 
2840 
2860 
2880 
2900 | 
2920 | 
2940 | 
2960 
2980 | 
3000 | 
3020 | 
3040 | 
3060 
3080 | 
3100 | 
3120 | 
3140 
3160 | ¢ 
3180 
3200 | 
3220 
3240 
3260 
3280 
3300 
3320 | 580 
3340 | 620 
3360 | 662 
3380 | 703 | .442 |: | 0.6) | 83. 
3400 | 744 | 433 : 21 S88 | SE 30.5 
J, a oem ne ae ee ee ee | 1100. .8 | 352. 168.3 
Calculated percentage transmission | 87 15.3 
Observed percentage transmission | 14.8 | 5. 83 
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The relative spectral emission from the lamp, for a given tempera- 
ture, is the product of the relative blackbody radiation intensities, 
the emissivity of the tungsten filament (each for the given tempera- 
ture), and the transmission of the fused-quartz lamp envelope. In 
table 2, the blackbody radiation and the emissivity values for a 
temperature of 2,737° K (which was the mean true temperature of 
the lamp filament in this case) are given, together with their product. 
The transmission of the quartz envelope has been neglected in this 
case, and for all practical purposes this product represents the relative 
spectral emission of the lamp for that temperature, since the trans- 
mission of the fused-quartz envelope varies but slightly throughout 
the spectral range in which this phototube is sensitive. 


IV. APPLICATION OF THE LAMP TO THE SPECTRAL 
CALIBRATION OF PHOTOTUBES 


The use of a continuous spectrum source in conjunction with filters 
to supplement the spectroradiometer in the determination of the 
spectral response of phototubes was discussed in detail in a previous 
publication [1]. Since the tungsten-in-quartz lamp is calibrated by 
a single measurement of the mean filament temperature, its spectral 
emission is readily determined and kept constant by controlling the 
voltage. As the result of heat conduction along the filament to the 
supporting elements within the lamp, and because of variations in 
filament thickness and absorption of radiation from other elements of 
the filament assembly, there is always a certain amount of variation 
in the temperature of different sections of the filament. Hence it 
becomes necessary to obtain a mean, or average, value for the filament 
temperature. This may be obtained by different methods. but for 
the purposes for which this lamp is to be employed it is probabl y 
most practically obtained by observing the color temperature * of 
the total visible radiation from the lamp and reducing this value to 
true temperature [6, 16]. 

Some calculations, in which the published values of emissivity were 
used [5], on filament sections having assumed differences in tempera- 
ture of the magnitude for the type of lamp described in this paper, 
indicate that mean true temperature values obtained from color- 
temperature measurements also accurately apply in the ultraviolet 
spectral range 2500 to 3500 A. 

The authors have considered the use of a ribbon instead of a round- 
wire filament in order to reduce the end and loop-cooling effects, but 
because of difficulties in keeping different sections of the filament in 
the same plane, it appears that the round wire is more desirable from 
the standpoint of the stability of the radiation received by the radi- 
ometer. The end cooling effects are reduced to a low value if the 
shields cover ¥% to % inch of the filament ends and loops. The tempera- 
ture variation between different sections of the filament could be 
reduced by a wider spacing of the filament loops. 

Since it is usually impossible to irradiate evenly the entire surface 
of the cathode of the phototube in the spectroradiometric method of 
calibration of the relative spectral response, the response curve ob- 
tained applies only to the portion of the phototube irradiated, and for 


‘These color temperature measurements were made by Deane B. Judd, of the National Bureau of 
Standards. 
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accurate work it becomes necessary to correct this curve so that jt 
will apply to the entire phototube surface as used in practice. 

Phototubes, as pointed out in the previous paper [1], often possess 
a long wavelength ‘tail’ in their relative spectral-response curves 
whic hi is not detected in the spectroradiometric calibration of the rela- 
tive spectral response, but which in the use of the tubes may produce 
very large errors, especially if the source of radiation has a rela tively 
high emission in this spectral region. Not only is the presence of 
long wavelength “tail” to the spectral- response curve of a phototube 
simply detected by means of this lamp in conjunction with a set of 
calibrated filters, but the magnitude of its effect may be closely evalu- 
ated, so that the phototube may be used in precision work. 

Since, as illustrated by the data of tables 3 and 4, the true photo- 
electric-response curve may, for all practical purposes, be much ties: 
ent from the spectroradiom etrically observed curve, especially nea 
the long wavelength limit of the tube, extra precautions in the use of 
double ~monochromators, etc., in the determination of phototube 
response curves would be useless effort unless the data were supple- 
mented by other means, for example, by observations using a tungsten 
lamp and calibrated filters. The use of a tungsten lamp and calibrated 
filters to supplement the results obtained appears to be required with 
even the best of spectroradiometers if high precision is to result. It 
would be very much better to use monochromator data, although it 
may not be of high accuracy, supplemented by lamp and filter measure- 
ments, than to rely solely on the best double monochromator. 


raBLE 3.—-Results of comparative tests of two phototubes for agreement beiween 
calculated and observed filter transmissions 


[The data for tube A were taken from table 2, and show close agreement between observed and calculate 
transmissions. The data for tube B show a large discrepancy, although the response curve was obt: sine 1 
with great care on a good speciroradiometer. ‘The calculations marked I were made upon the basis 
smooth spectral-response curve drawn through the spectroradiometrically observed points. Calculations 
II were made upon an adjusted curve (see table 4). Tube B appeared to have a slight amount of “‘/atigue, 
which ——w accounts for part of the discrepancy between the calculated and observed filter trans- 
missions 











Photosmbe4 (20 | Phototube B (no filter) | Phototube 2 (Nite Ate 


| r a 
Cale.I | Obs. | Cale.I | Cale. II | Obs. | Cale.I | Cale. II 


Filter 





75. 1 


53.8 


57.9 | 38. 2 | 4. 34.0 | 75.6| 63.5 
32. 0 19. 3.7 | % 53.8 | 31.4 


15.6 


15.3 | 8.8 ‘2 | 61 313 ro 34.8 


_— 3.3 | .1 | : 13. 4 | 1.0 
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TABLE 4.—Adjustment of phototube spectral response 


(The values indicate the adjustment in the relative spectral response of phototube B required to give filter 
transmissions which approximate the observed values when using the tungsten-in-quartz lamp and 
calibrated filters, see tables 2 and 3] 


| Spectroradio- 

| metrically Adjusted 

| Wavelength observed response, 
response, curve II 
curve I 


A 


2600 
2640 
2680 
2720 
2760 
2800 
2840 
2880 
2920 
2960 
3000 
3040 
3080 
3120 
3160 
3200 
3240 
3280 
3320 
3360 
3400 


Cowan 


Os eho com Ons 
o 





ee ere a ae . ee 





By measuring the integrated transmissions of a set of calibrated 
. t) D 


filters with this lamp as a source, and by calculating the integrated 
transmissions of the same filters from a knowledge of their spectral 
transmission, the spectral-energy emission of the lamp, and the spectral 
response of the phototube (based upon an arbitrary curve drawn 
through the values for the wavelengths at which the response was 
determined spectroradiometrically), the shape of the phototube 
response may be altered until the observed and calculated transmis- 
sions are in agreement. 

For all practical purposes, the number of filters can be increased by 

n—1) by using combinations of each filter with the one having the 
nearest long (or short) wavelength cutoff. In certain cases more 
remote combinations can also be employed to advantage. This is 
especially true where filters having bands of selective absorption are 
emploved in evaluating radiation over a long wavelength range of the 
spectrum. The particular problem in hand, together with the 
characteristics of the available filters, will determine the most useful 
combinations. A great amount of labor may be saved by using a 
combination filter rather than calibrating a new one where the com- 
bination gives the desired spectral transmission. 

In this method of determining the relative spectral response of a 
phototube the entire electrode of the cell is fully and evenly irradiated, 
so that any local variations in the relative spectral response are 
averaged into the final response curve for the cell as used in practice. 
This is important, since in the manufacture of cells a uniform cathode 
surface is rarely obtained. Differences in the thickness of the sensi- 
tive layer or the presence of impurities modify the relative spectral 
response of one area of the receiving surface relative to another ad- 
joining area in the same phototube. 


524586—43——3 
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Having obtained a careful evaluation of the spectral photoelectric 
response and the transmissions of the available filters, the problem of 
the measurement of ultraviolet radiation from a given source requires 
a certain amount of information concerning the spectral quality of 
the source under investigation. If the source is a temperature 
radiator, as, for example, an incandescent-filament lamp of unknown 
filament temperature, a smooth curve may be assumed and adjusted 
to its true relative spectral value through a set of observations and 
calculations on the filter transmissions. For a discontinuous source. 
as, for example, a metallic arc, the wavelength and a knowledge of 
the relative intensity of the emission lines are necessary. In the case 
of a source of intermediate type, as a carbon arc or the sun, where 
there are emission bands with intervening regions of low intensity, a 
knowledge of the shape of the emission curve is essential. 

In the calculations on the sun [1, 8, 9] a smooth curve was assumed 
(because the true shape of the curve was unknown), and tbe results 
are given in terms of the equivalent radiation for a definite waye- 
length interval on the basis of that assumption. This has a short- 
coming, in that it is impossible to evaluate accurately the intensity of 
the radiation in any definite wavelength interval, since ultraviolet 
solar radiation is not distributed according to a smooth curve but 
is concentrated in definite bands with relatively low intensities at 
intermediate wavelengths. The accuracy of any evaluation of this 
type will therefore depend upon the wavelength interval chosen and 
upon the range of spectral sensitivity of the phototube employed. 

Therefore, it at once becomes obvious that to evaluate ultraviolet 
solar radiation accurately in a definite wavelength interval an approx- 
imate shape of the true spectral solar radiation emission curve should 
be obtained. Since it is impossible to obtain this curve by means of 
an ordinary spectroradiometer (because of scattering within the 
instrument), it would be a valuable contribution in this field for some- 
one having available a quartz lens and prism double monochromator 
to make this determination for (different air masses) several altitudes 
of the sun. Photoelectric filter radiometry in solar radiation studies 
will, with this accomplishment, be put on a practical scientific basis. 


V. APPLICATION OF THE LAMP FOR THE TOTAL 
RADIATION CALIBRATION OF PHOTOTUBES 


A practical method of using a standard quartz mercury-arc lamp for 
calibrating the total response of phototubes, together with the asso- 
ciated amplifier or other indicating or recording instrument, has been 
described in detail elsewhere [8, 9, 10]. The herein-described tung- 
sten-filament-in-quartz lamp may be used in a similar manner and 
may have advantages in some cases over the mercury lamp. The 
radiation intensity is more constant, both in relative spectral quality 
and in total intensity, and its operation requires less care in the con- 
trol of room temperature and circulating air currents. Furthermore, 
since the tungsten lamp emits radiation throughout the range of 
sensitivity of the phototube, the calibration applies to the total 
integrated curve rather than to a few isolated wavelengths. If the 
cell is employed in work involving a more or less continuous spec- 
trum, or one which is at least different from the mercury arc, the 
advantage of using this type of spectrum in the calibration 1s obvious. 
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Since the lamp envelope absorbs the long infrared radiation emitted 
y the filament and emits some infrared radiation, of very long wave- 
length, characteristic of the temperature of the quartz envelope, the 
lamp is more easily calibrated, for absolute radiation intensity per 
unit wavelength at a fixed position relative to the lamp, by employ- 
ing a 1-centimeter water cell over the thermopile. Phototubes are 
insensitive to this infrared radiation, and a water cell is not required 
in the use of the lamp with them. 

If the sensitivity of the phototube is confined to a narrow range of 
the spectrum, as in the case of the two titanium tubes referred to 
ibove (see tables 2 and 3), additional filters may be employed to 
evaluate the limited range of the spectrum of the lamp radiation cover- 
ing that to which a particular phototube is sensitive. Since the 
relative emissivity of tungsten is known none too well over the total 
spectral range from 2300 to 15 000 A, it is obviously an advantage to 
limit the range of calibration of the lamp to the spectral range over 
which the phototube responds. 


1 
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EFFECT OF TURBULENCE AND CHANNEL SLOPE ON 
TRANSLATION WAVES 


By Garbis H. Keulegan and George W. Patterson 


ABSTRACT 


This paper is the second of a series dealing with the motion of flood waves and 
other waves of translation in open channels. The first paper considered waves 
controlled solely by inertia forces; the present one is an analysis of the combined 
effects of turbulent friction and inertia. The basic equation ‘of motion for gradu- 
ally varied unsteady flow in prismatic channels is derived from fundamental 
ans iciples. The effect of the velocity distribution in the original undisturbed 

urrent on the motion of short waves is investigated, and the effects of wave 
height, curvature of profile, and fluid friction on the celerity of a wave-volume 
ele ment is analyzed in detail. The deformation of a straight sloping front and the 

iange of height of an abrupt wave front is treated. Special emphasis is laid on 
te sce of negligible curvature and practical methods of handling engineering 
problems arising in connection with the operation of locks or hydroelectric canals 
are given. 
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LIST OF SYMBOLS 


Symbols used in only one section are omitted unless they are of general signifi- 
cance. Important quantities in hydraulic equations are indicated as functions 
All equations are valid in any consistent set of units. The significance 


of z and t. 


of barred variables is explained in section II-1. 


A(z, t) 
b 


exp ( 
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fi 
fa 


g 
H(z,t) 
Ho 


h(a, 


t) 

1 
ky 
ks 


l, m, n 
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Pa 


Pzz + 


Res, 


q 


R 
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cross-sectional area of channel. 

width of free surface. 

exponential function. 

base of natural logarithms. 

arbitrary functions, determined by boundary conditions. 

substitution variable, (3w9—5Up)gi/3 Uo(wo— Us). 

substitution variable, o/4 or — do/4. 

acceleration of gravity. 

depth of water. 

undisturbed depth of water. 

increase of depth due to disturbance. 

slope of channel bed. 

substitution variable, 1—3Up/wo. 

substitution variable, 1— U /w». 

direction cosines. 

pressure. 

atmospheric pressure. 

generalized stress components. 

discharge per unit width of channel. 

turbulent stress components. 

hydraulic radius. 

viscous stress components. 

time. 

mean velocity in channel cross section. 

undisturbed mean velocity. 

increase in mean velocity due to disturbance (sections III-3, ff.). 

x-component of velocity (mostly before section III-3). 

velocity fluctuation due to turbulence. 

y-component of velocity. 

velocity fluctuation due to turbulence. 

z-component of velocity. 

velocity fluctuation due to turbulence. 

Cartesian coordinate, specifically, parallel to channel in bottom 
plane. 

Cartesian coordinate, specifically, lateral to channel in a hori- 
zontal plane. 

Cartesian coordinate, specifically, drawn upward with origin at 
bottom of channel. 

Boussinesq coefficient of velocity distribution. 

eddy viscosity. 

coefficient of channel friction. 

coefficient of channel friction, H= Hp. 

viscosity. 

density. 

partial volume of wave per unit width. 

partial volume of wave per unit width. 

mean shear per unit area on channel wall. 

arbitrary function. 

wetted perimeter. 

arbitrary functions. 

velocity of a wave-volume element. 

velocity of a short wave, moving without change of form. 


I. INTRODUCTION 


This paper is the second of a series dealing with the motion of flood 
waves and other waves of translation in open channels. The purpose 
of the series is explained in the introduction to the first paper [1].’ 


1 Figures in brackets indicate the literature references at the end of this paper. 
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\fany of the significant properties of translation waves may be ob 
tained by assuming that the motion of fluid particles due to the wave 
disturbance is irrotational [1]. Experimental observation of these 
properties indicates that they are in good agreement with theoretical 
predictions, provided that the waves are moving in still water. If, 
nn the other hand, the disturbances are produced in a liquid which 
is initially in motion, as in actual currents in an open channel, some 
disparity will be seen between observation and theory. Bazin’s ex 
periments on the velocity of propagation of a solitary wave in still 
water or in a running current demonstrate these points [2]. 

The divergence between observation and theory, when the wave 

s moving in a current, can be explained in two ways. First, irrota- 
tional or potential theory supposes that the original undisturbed 
motion is such that all the particles in a cross section move with the 
sue constant velocity, Up, which may of course be zero. This 
condition is not fulfilled for flow in actual channels, and it thus be- 
comes necessary to investigate the effect of the distribution of velocity 
in a cross section on the velocity of propagation of waves. An analysis 
of this effect has been carried through in great detail by Boussinesgq [3]. 
\With a view to simplifying the presentation and the solution of the 
problem, the authors give a new derivation (section III-5). Al- 
though the results arrived at here are different in form from those of 
Boussinesq, it is believed that the difference has a qualitative rather 
than quantitative importance. 
This difference is due to two distinct forms of the equation of 
motion of unsteady flow in a channel, one adopted by Boussinesq and 
the other by the authors. The present method applies the theory of 
Reynolds stresses to the problem of the unsteady mean flow in a 
channel, without assuming any specific relation between the stresses 
and the velocity rates existing at points within the channel cross 
section. The effect of these stresses is evaluated in terms of the total 
friction at the wall, which is the important concept for mean flow in 
open channels and is usually available from observations on uniform 
flow. Boussinesq’s method on the other hand is an extension of his 
own theory of turbulent flow, and specific relations between the 
turbulent stresses and the space-rates of the local velocities are as- 
sumed and introduced into the derivations. As the concept of 
Reynolds stresses is the starting point of the present analysis, the 
theory of these stresses is given in considerable detail (section II-1). 
The equation of motion of a perfect fluid is the point of departure in 
interpreting the relation between the stresses and velocity fluctuations. 
This novel procedure simplifies the presentation considerably. It is 
hoped that the equation of motion for gradually varied unsteady flow 
in prismatic channels derived in this paper, eq 47, is sufficiently 
rigorous to merit the attention of river engineers. 

In the second place, the irrotational theory of translation waves 
ignores the effect of turbulence, which is one of the forces controlling 
the motion of the waves. Consider a short intumescence moving in a 
current. In the initial stages of motion the behavior of the wave is 
not affected by turbulence, since the wave is assumed to be short. 
When the wave has traveled a considerable distance the cumulative 
effect of friction becomes appreciable, particularly at the rear of the 
wave. In fact, experience shows that the front of the wave will 
travel with a velocity characteristic of irrotational waves, whereas at 
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the other end of the wave, the volume elements will move with a 
smaller velocity, the exact law of which is not as yet clearly understood. 
The wave thus flattens out considerably at the posterior end, with a 
consequent elongation. The tendency of short waves to deform in this 
manner is clearly demonstrated in Horton’s experiments [4]. 

When the motion of the waves is controlled principally by friction, 
as in the quasipermanent regimes of rivers, the corresponding theory 
of translation waves is simple. But when it is supposed that the fric- 
tion and inertia effects are of like magnitude, analysis is confronted 
with a difficult task. The change of shape of a short intumescence 
during the entire journey of the wave in a current would be accurately 
described if a complete solution were available, which is not the case. 
However, some progress in this direction is possible if it is assumed that 
the effect of turbulence is small, requiring only a secondary correction. 
Such a theory has been worked out by Boussinesq, who assumed that 
the velocity of the wave-volume elements differs little from the theo- 
retical velocity of waves without friction [3]. The main purpose of this 
paper is to consider this theory and the corresponding solutions for 
waves of negligible curvature. In giving Boussinesq’s solution, slight 
changes have been made, since it was desired to adopt Manning’ s law 
of friction in open channels [5]. Favre has employed the Boussinesq 
integrals to describe the motion of waves occurring during the opera- 
tion of locks [6]. Because of the engineering significance of the problem 
we have also reproduced these applications (section IV-9 and IV-12). 

The applications of the theory are presented in sections III and IV, 
which have been entitled first-order and second-order theories of 
wave propagation. The first-order theory neglects the vertical acceler- 
ation of the water particles and assumes that the combined effect of 
friction and velocity issmall. It is found that this assumption requires 
that the waves be short. The second-order theory considers the effects 
of vertical acceleration but neglects the variations of the horizontal 
velocity components. As in Boussinesq’s investigations, the concept 
of the velocity of a wave-volume element plays a very basic role. 
Equation 116 is the general expression for this velocity. Using this 
same formula, the authors have been able to study the gradual de- 
formation of a straight sloping front and also to compute the reduction 
of height of an abrupt wave front. 


II. FUNDAMENTAL PRINCIPLES 
1. REYNOLDS EQUATIONS OF MEAN MOTION 
In the Reynolds theory of turbulence [7; 8, p. 638; 9, p. 364; 10 
p. 191] the velocity components u, v, w, in the directions z, y, 2, 
een are regarded as fluctuating functions of the variables 
x,y, 2,t. The fluctuations are frequent and rapid, and the periodicity 
of the changes is irregular. Denoting the mean values of wu, v, w by 


u, v, w, and the mean value of the pressure, p, by p, we hae for the 
instantaneous values 


ms , 

u=utw’, 
v 

=v+v ' 


w=w+w’, 





p=pt+p’, 





Translation Waves in Open Channels 465 


where wu’, v’, w’, and p’ are the fluctuations whose mean values u’, v’, 


w’, p’ are equal to zero. The mean values may be temporal; that j is, 
taken at a given point (x, y, 2) over an interval of time of duration 


At; thus, for example, 
— 1 ite 


Or again the mean values may be iii that is, taken at a given 
instant ¢ over a space of volume AV; thus, for example 


- 1 : 
yx ap for udV, (3) 


For the present it will be supposed that the averages are temporal. 
A difficulty arises when taking temporal means in unsteady flow. In 
such cases it is necessary to make the time interval At long enough to 
include a great number of fluctuations, but short enough to exclude 


all but negligible variations of u, v, w, p. 
There are certain rules for taking the averages. Expressed sym- 
bolically [9, p. 366] 


and 





where a and 6 are functions of z, y, 2, t; s represents any one of the 
independent variables, and a bar over a quantity represents an average 
of the kind just described. These rules are based upon assumptions 
that are not exactly applicable to the flow of liquids. The errors, 
however, are negligible when the fluctuations within the time interval 
At are sufficiently numerous. 

The specific contribution of Reynolds’ analysis is the idea of in- 
terpreting the dynamic effects of the fluctuations w’, v’, w’, as apparent 
tractional forces or stresses. The nature and the magnitude of these 
stresses can be readily obtained starting with the equations of motions 
of a perfect fluid. The dynamic equations in question are 


du p sf 
P'dt 32 toX, 
dv___ op » 


dw_ _Op 
p dt 52 1 PS: 





where the operator d/dt, in ene of the local or is 


ad. «0 
dt =S +g tes, bw :. 
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Together with these, there is the equation of continuity, assuming 
that the liquid is incompressible, : 
Ou , OV , Ow __ 0 me 
Ox Tay te ow 4) 
The liquid will always be assumed to be incompressible, in this paper, 
In eq 5, p is the pressure, p the density of liquid, and X, Y, Z the 
components of the force of gravity per unit mass of liquid in the 
direction of z, y, 2. These relations are valid at any instant of time 
independent of the fact that the quantities, u, v, w, and p are fluc- 
tuating functions of time, provided that the velocities are not dis- 
continuous, and if discontinuities exist, they occur only occasionally 
during the time interval At. . 
Replacing u, v, w in eq 7 by their respective values from eq 1, 
averaging each term, and using the last three relations in eq 4, together 
with the fact that u’=v’=w’=0, we obtain 


Ou ,Ov , OW 
dr dy Tez” 8) 


which states that the dilatation of the mean values of the velocity 
components vanish. It then follows from eq 1, 7, 8 that 


ou’ 
or rs (9) 


The equations of motion, eq 5, are treated in a similar manner. 
We use 6, substitute for u, v, w their values from eq 1, multiply eq 9 
by pu’, pv’, pw’, respectively, and add to the left-hand side of each 
equation. Each term is then averaged in accordance with eq 4, 
assuming p to be constant, that is, that the fluid is homogeneous, 
obtaining - Me 

du_ Op — = 
Odi ant OATS, TS TOs” 


dv_ Op, Rr Ry Ry 
Fat oy Ft T Se ey Tt Se’ 


0 __OP 7 Pas Ry, Pe 
Pat ae MTS T3y oz” 





where, in this case, 
d 0,-0 ,-d0 .-0 
dt ott “dat Oy t M9, 
and 
R,.=—pu'u’; R,, =R,,=—pu'v’, 
Ry=—pv'’; R,=Ryz=—pu'w’, (12) 
Sompeey — 
R,=—pw'w’ ;R,,=R,=—pv'w’. 


We may refer to eq 10 as the Reynolds equations of motion for 
turbulent flow. 
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To interpret the above results we write the general equations of 
motion of a deformable continuum under the action of the external 
body forees X, Y, Z and surface tractions of an unspecified physical 

origin: 
du __ X4 Be 4s 
‘a 
dv 
‘a 


Oz 
ie Oz’ | 
OPyy OP ev, 
Oy Oz’ 


oy 
pY Ts Pes 


dw 


pag mee ee $e Bee 4s 





where the operator d/dt is the one Pre in eq 6. 

It should be noted that these equations are merely Newton’s law of 
motion as applied to any deformable continuum, and thus they apply 
to a wide range of phenomena. The left-hand members represent 

components of the rate of change of momentum per unit volume of the 

medium, whereas the right- hand members represent the components 
of the forces acting per unit volume. The nature of the medium 
determines the functional relation between the nine stress components, 
Per, Pyty » » . . . , and the other variables occurring in the equations 
and the intrinsic properties of the medium. Hydraulic equations 
established by integrating these equations over a cross section of the 
fluid normal to the flow are frequently referred to as momentum 
equations. 

The nine stress components completely determine the forces acting 
on a small surface dS at the point (z, y, z). If a normal is drawn to dS 
away from the bedy of which dS forms a part of the surface, and if /, 
m, n are the direction cosines of the normal to dS, then the z-, y-, 2- 
components of the surface tractions on the body at the point (z, y, z) are 


F,=lpyet+mMpyzt+Npizy 
Fy=lpryt+Mpyt+npx, (14) 
ant 
F,=lp..+ MP y2+ NP zz- 


Comparing eq 10 and 13, we see that 
Pr = —p+Rrs, Pay=Pys= RyRy; 
Py=—pPtR,,, Pu=Pu=Rh,=R; 


P22 —p+R.,, Py=Py=R,=R,;; 


and thus the turbulent mean flow may be ee as ordinary stream- 
line flow where the stresses Prz, Pry... . , are identified with the 
stresses R,,, R,, ..... These apparent Noi of turbulence are 
called Reynolds stresses, and represent mathematically the transport 
ie momentum across a surface caused by the velocity fluctuations 
10, p. 192]. 

If the flow instead of being turbulent is viscous in character, the 
stresses are of molecular origin. The stresses Pzz, Pry, - . - are ‘then 





468 Journal of Research of the National Bureau of Standards 


identified with the viscous stresses S,z, S,,, ... in the following 
manner, 


Paes PA Ser, Pry = Piz = Sry? = Syz} 
Pu=—Prt Sys, Pre Pez= =S2,=S8::; | 
Pes= —P+Wr15; Py2=Pay=Sye= Sry. 


If the liquid is incompressible, these viscous stresses become 
Ov , Ou 
S —_ S, — = 
ry Mu Br toy)’ 


se . ow =) 
=Sr=nl 9, F895 }} 





i ag 
S,.=S,= W\ Set dy 


where yu is the viscosity of the liquid. 

It is perfectly feasible to carry through the derivation of the Reyn- 
olds equations of motion, eq 10, using the viscous-flow interpretation 
of eq 13 as a starting point, and finally neglecting the viscous stresses 
in comparison with the Reynolds stresses, which is justifiable on the 
basis of experimental evidence [11, p. 186]. We have not adopted 
this method, since the results are identical with those obtained here 
by initially assuming a perfect liquid, that is, by neglecting the vis- 
cous stresses at the beginning rather than at the end, and since it is 
more complicated than the me ~. used here. 

Later, in sections III-1 and IV-1, we shall derive hydraulic equa- 
tions for mean flow in open cae using eq 10 as a starting point. 
For this purpose, it is sufficient merely to suppose that P,,, Z. e 
exist, there being no necessity for determining them as functions of the 
mean velocity gradients, because when the Reynolds equations of 
motion are integrated over the cross section of the channel, the 
quantities R,,, R,,, . .. either disappear or are expressed in terms 
of the wall friction, 7», which is defined in connection with eq 45. 
Boussinesq, however, expressed R,,, R,,, . . . in terms of the velocity 
gradients. He simulated eq17, merely replacing S,;, S.,, ... by 
R,,, Ry, ... ju, v, w by u,v, w; and wp by e. The quantity ¢ was 

called the coefficient of internal friction by Boussinesq [3, p. 46]; 
present usage seems to favor the term eddy viscosity, since ¢ has the 
same dimensions as u [11, p. 185]. 

When the channel is rectangular, Boussinesq assumed that « is 
constant for every point in a given cross section, but in a channel of 
constant roughness, directly proportional to the product of the wall 
velocity and the depth of the liquid. Such assumptions form the 
foundation for the developments in his ‘‘Essay.”’ The validity of his 
theory is thus open to question, since the assumptions concerning « 
are widely at variance with present experimental knowledge. 

Having established the basic partial differential equations governing 
turbulent flow in open channels, we now turn our attention to the 
boundary conditions. Consider first the condition to be satisfied at a 
fixed solid wall, which we shall suppose is covered with uniform 
asperities of height k. A pillbox is erected on the surface, with height 
h, just sufficient to clear the asperities, and radius r, large with respect 
to 2h. (See fig. 1.) 
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Figure 1.—A pillbox enclosing the asperities on a channel wall. 


Since the liquid is incompressible, the net rate of flow across the 
surface of the pillbox must be zero. We neglect the discharge from 
the cylindrical side of the pillbox, since its area is small compared with 
that of the circular end. There must be at least one point on the end 
of the pillbox for which 

lut+-mov+nw=0, 


where 1, m, n are the direction cosines of the outward-drawn normal. 
Averaging the terms in this equation and assuming that wu, v, w vary 
slowly in an interval of length 7, we have for all points on the fixed 
surface 

lut+meo+nw=0, (18) 


which is the desired boundary condition for a fixed solid wall. If the 
wall is smooth, the argument is repeated, replacing the height of the 
asperities by the thickness of the laminar boundary layer. 

Next we consider a boundary condition which applies to the free 
surface. Let F (2,y,2,t)=0 be the equation of the free surface. As 
before in the case of the velocity components and the pressure (eq 1), 
the surface is considered to fluctuate rapidly and irregularly about 
a mean value, and we write 


F+F’, (19) 


where F is the mean value of F, and F’ is the fluctuation whose tem- 


poral mean value #” is zero. 
Since particles on the surface remain there, we have the condition 


OF, QF DF OF 
yy Oat ay te oz” 


where u,, V;, Ws, are the velocity components at the free surface. If 
we replace u,, vs, ws by expressions obtained from eq 1 and F by its 
value from eq 19 and average the terms, the result is 


-oF = 7 
eee ioe ape “=0 (20) 


provided we neglect the terms w’,(0F’/dz), .... These products 
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may be ignored, since either the correlation between w’,, OF’ ou, 
etc., is zero or the term is small in comparison with those retained, 
On the basis of this assumption, the above equation is a boundary 
condition obtaining at the free surface of turbulent flow. 

In order to simplify the notation from this point on, 4, v, w, p, 


etc., will be written for u, v, w, p, etc. Any further use of the bar to 
indicate an av erage will be in the macroscopic, rather than th; 
microscopic, sense. 


2. EQUATION OF CONTINUITY FOR MEAN FLOW IN CHANNELS 


Up until this point we have derived equations which are valid for 
turbulent flow in general. We now derive an equation which applies 
to any tube of flow with a cross section which varies both in space and 
in time. This tube is assumed not to close on itself, and since it is a 
mathematical representation of a flowing channel, it is drawn as such 
in figure 2. The cross-sectional area of the channel, normal to 2, will 
be represented by A=A(z,t), the average value of wu over the cross 
section by U=U(z,t). 














FiGURE 2.—Surface configurations at two different times in an open channel. 


The positive sense of x is chosen so as to make U positive. The 
actual direction of z is arbitrary, but for convenience it should approxi- 
mate the direction of the channel. 

The flow into the stretch Az during the time At is 


to+At 
f. U (aa,t) A (2a,t) dt, 


while the flow out of the same stretch is 


to+At 
j, U (a+ Az, t) A (2o+-Az,t)dt. 





W 


—imakh «ab abit A. 
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Furthermore, the volume of liquid contained within the stretch at 
time fy 1s 


TotAr 
A(z,t))dz, 


Zo 


while the volume in the stretch at time t+ At is 
Zo+Ar 
A(z,to+ At)dz. 
Zo 
Since the fluid is incompressible, net outflow plus increase of volume 


equals zero, or 


ot At 
"aot da, f)A (a+ Aa,,t) — U yt) A (rat) 


ai 
e 


To+Ar 
nm j, [A(2x,to-+ At) —A(z,t,)]dz=0. 


If we apply the theorem of the mean for integrals, divide by AzAt 
and then let Ar and At approach zero, we obtain 


Orr 4) 0A _ (21) 
5 (VA) +3 = 


which is the equation of continuity for mean flow in channels. Since 
no assumptions have been made regarding the shape of the channel, 
this equation applies to open channels in all generality. 


III. FIRST-ORDER THEORY OF WAVE PROPAGATION 


1. THE EQUATION OF MEAN FLOW IN PRISMATIC CHANNELS 


Consider a prismatic channel and select as the origin of a right- 
handed system of Cartesian coordinates tbe lowest point of any con- 
venient cross section (see fig. 3). Take the z-axis as the line parallel 
to the channel (not necessarily horizontal), with positive sense in the 




















Figure 3.—System of coordinates in a prismatic channel. 


direction of flow, the y-axis as the horizontal line, and the z-axis 
positive upward (not necessarily vertical). The channel slope, 2, is 
the tangent of the angle 6 from the positive direction of the z-axis to 
the horizontal, measured in the counterclockwise direction, viewed in 
the direction of increasing y; in other words, a positive slope is down- 
hill. These conventions will be adhered to throughout the paper. 





472 Journal of Research of the National Bureau of Standards 


We assume that the angle of slope of the channel is so small that 
sin @ is sensibly equal to tan @ and cos 6 to unity, hence X and Z in 
eq 10 are gi and —g, respectively. 

Again, in order to simplify the three equations of motion, eq 10, 
we make two distinct sets of assumptions, the first involving the 
terms on the left-hand side and the second involving the terms 9n the 
right-hand side. First, we assume that the second derivatives of the 
cross section, A, and of the mean velocity, U, and the squares and 
products of their first derivatives, may be neglected in comparison 
with the first derivatives. Such a flow is referred to by Boussinesgq 
as gradually varied [3, p. 242]. In general, this requires that the 
sides of the channel be not much flatter than a 1:1 slope, that the slope 
of the wave profile be small compared to unity, and that the curvature 
of the wave profile be small compared to (1/U)(0U/0z) or (1/U%) 
(OU /ot). Asa result of these restrictions, dv/dt and dw/dt are small 
with respect to du/dt, and the left-hand sides of the last two of eq 
10 vanish. 

Secondly, we make certain assumptions concerning the Reynolds 
stresses: (1) u’u’=v'v’=w"'w" hence R,,-=Ry=R.2; (2) since u is the 
predominant velocity component, the derivative of the correlation 
between the fluctuations v’, w’ may be neglected in comparison with 
the derivatives of the correlations between wu’, v’, and uw’, w’, which is 
to say, that the derivatives of R,, may be neglected; (3) since the flow 
is gradually varied, the z-derivatives of R,, and R,, may be neglected. 

On the basis of these two sets of assumptions, then, eq 10 reduce to 


Ou, Ou, Ou, OU . te 10 10 
— me —_ —— — f j— a | —-- —- 99) 
ya "he I vt u ~~ gj > Sr? Mis tant, 5 alten (22 


r) 
0=5,0 = 


0=+ a oN (24) 


The simultaneous solution of eq 23 and 24, with the dynamical 
boundary condition p=p, at z=H, is 


P—Ri2=Patpg(H—z), 


where H, a function of z and ¢ only, is the z-coordinate of the free 
surface, and p, is the atmospheric pressure. These symbols will be 
consistently used to refer to these quantities. According to the above, 
in gradually varied flow the pressure is hydrostatic everywhere. 
Differentiating the above equation, it is found that 

O(p Rex) _ gH 


Or J Or 


Each term in eq 22 is to be integrated over the cross section A at z 
In this way, the resulting equation of motion will involve only the 
mean flow, U, the cross-sectional area, A, and certain other quantities 
to be given later. We thus eliminate the variables y and z from con- 
sideration, since all of the quantities occurring in the final result will 
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be functions of x and ¢ only. Such an equation is called a hydraulic 
equation, and bears the same relation to circuit equations in electricity 
as the general equations 10 do to the electromagnetic field equations. 

Let us first evaluate the integrals arising from the terms in the left- 
hand side. The mean flow, U, in a cross section, is implicitly defined 


by 
UA | wl’ 
A 


We write the above as a double integral 


, H v2 
UA- f, az{ udy, 
0 Jn 


where y2 and —y, are the y-coordinates of the two points on the wetted 
perimeter having the common distance z from the z, y-plane (see 
fig. 4). Obviously, y. and —y, are functions of z only, for a given 
prismatic channel configuration. Differentiating both sides of the 
above equation with respect to t, we obtain 


O(U"A) __ OU 44 oH (” i 
ue iat” Oi 


where 62 and —6, are the y-coordinates of the two points on the wetted 
perimeter at the distance H from the a, y-plane. 


Z 


A 









































Yry 


Figure 4.— Double integration with respect to y and z; cross section, looking upstream’ 
524586—43-—_-4 
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We now define u, to be the mean value of the surface velocities. 
us, over the width b=b,+ 5, of the free surface: 


a bs 
ub= f u,dy, 
—h 


and the next to the last equation becomes 

rr Ou,, ,—OA 

all A) = a at4At Us Ot ’ 
since b(0H/dt)=O0A/ot. Introducing the velocity coefficient a;, implic- 
itly defined by 
a,U=u,— U, (29) 
eq 28 may be written 


ou ou OA 
st A=A > U 7 


We now turn our Cees to the second term on the left of eq 22. 
We define another coefficient of velocity distribution a, as the ratio 
of the mean of the square of velocity to the square of its mean, both 
means being taken over the cross section, thus 


ara | wdA. (31) 


This coefficient is the only one which appears in the final equation, 
and was originally introduced by Boussinesq. It also plays an 
important role in the equation of steady flow in open channels [12] 
and must not be confused with the so-called Coriolis coefficient, which 
is frequently denoted by a. Assuming that a does not depend on z, 
and differentiating both sides with respect to z, in the same manner 
as was done in oe eq 28, . have 


™ d(u® 
a2 (L “y= f 2% gA+ee 


Defining u,” to be the mean value of mi squares of the surface veloc- 
ities over the width, 6, of the free surface, 

— ba 

usb= udy, (32 


—b; 


then the previous equation becomes 
Drains 0(u”) — 0A 
hae = = = ee 
as (l 1) { - dA-+u, Sr’ 


since b(0H/Ox)=O0A/Ozx. Introducing a third coefficient of velocity 
distribution a2, defined by 
a2 UV?=u,?—aV?, 
we obtain from the nays two relations, 
ou 0A 


ou dA=al ‘AS —a,l* — 


(30) 


a : 
 uidy, 


Next consider the integrals of the third and fourth terms of the 
left-hand side of eq 22. For the third term we write 


‘yiA= { ie is sith 
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where the notation is the same as that used in the above two integrals 
» m . . . . 4 
see fig. 4). The integral with respect to y is integrated by parts, 
obtaiming 


ou sal ov 
[eSuaa at 3 [ (Up) 2— (Upro) ]dz — f aseaa, (35) 


where % and Up are the velocity components at the wetted perimeter 
and the subscripts 1 and 2 indicate the two sides of the channel, that 
is, the points (—y;,2) and (y2,z), respectively. 
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FicureE 5.—Double integration with respect to z and y; cross section, looking upstream. 


Again, for the fourth term 


ou bs H Ou 
fi» 5244 = fia wed, 


where 2 is the z-coordinate of a point on the wetted perimeter at the 
distance y from the 2,z-plane (see fig. 5). This time the integral 
with respect to 2 is integrated by parts, obtaining 


ou ba ow 
J Py * dA= i (WsUs— Won) dy — J us, 4A, (36) 


where Wo, Up are now evaluated at the point (y,2:). 

If dx denotes an elementary arc of the wetted perimeter, x, measured 
in the counterclockwise direction, and the direction cosines of the 
inward-drawn normal are 0, m, n (figs. 4 and 5) then dy=ndx and 
dz=—mdx, and{eq 35 and 36 become 
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. ov 
Unymdx— | u—dA, 
J diel JA Oy 


. Ou *be <7] Ow 
| w. dA | wu dy | wyuond x - ux—dA. (38 
Ja Of . = a 1 OZ 


We add these two equations together, noting that as a result of eq 18 
the integrals with respect to x vanish, and combine the integrals in the 
extreme —_ by means of eq 8 and 34 obtaining 


? " Ou ba ou OA 


Furthermore, since the equation of the free surface is 
F(2,y,2,t) =2z—H(z,t) =0, 


and since 0A/dt=b(0H/dt) anddA/or=6(0H/dz), the surface boundary 
condition, eq 20, becomes 

OA, OA 

bw,=— tu, 

ot Or 
Replacing w, in eq 39 by its value from the above equation, recalling 
that A is not a function of y, evaluating the integrals by means of 
eq 26 and 32 and eliminating 0U/dxz by means of eq 21, we have 


+f Te =(u,—aU pt +t U— al 


Introducing a final coefficient of velocity distribution, 
Ua,;=u,—alU, 


the previous equation becomes, using eq 33, 


\ 
* Ou : vou 7 0A OA ( 

1A 9 : . 
Ke >y' + Z Laz tel or f 
( 
We have now succeeded in evaluating all the terms on the left t 

of eq 22 in terms of functions of x and ¢ only, and proceed to the 

terms on the right. In view of eq 25, the first two terms become f 
merely ! 
t 
19 >H 

ans c 7) 
hata adiaa® Re) dA=nA- gA5, 
; ; ; t 
since // is not a function of y or 2 : 
The integral of the last two terms in eq 22 is transformed by Green’s 


theorem into a line integral around the periphery of the cross section: 


iffo, ,d vm 
“LSet Re p=} Gf Bydz—Redy. 
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Expressing the line integral in terms of x as was done for eq 37 and 38, 


4 O | 
= = t5oRee dA Pink, + nRJdx. 


Comparison with eq 14 and 15 shows that the integrand of the line 
integral is the z-component of the force which the fluid exerts on the 
channel walls and the free surface. If we denote this shearing force 
by 7, and assume that the shear on the free surface is zero, and define 
7, to be the mean shear of the fluid on the channel walls and bottom 


1 rad 
To - rd, 
Xe x 


fre .. 
LL SP +3 ie . |e 


Our choice of axes makes 7» positive, since 7) has the same direction as 
1’. We now define » the channel coefficient of friction, 


then we have 


Uj? 
T sell 


FI 


~ 


and eq 44 becomes 


tee . eo . AL *y 


Collecting the results of eq 30, 34, 41, 42, 46, and dividing them by 
A, the desired integral of eq 22 is 


pt ol’ OU »)UoA a " 
rt pratt pt I~) 4 op 9 OR” nied 


where R=A/x is the hydraulic radius of the channel. This is the 
desired hydraulic equation of motion for gradually varied unsteady 
flow in prismatic channels. It will be noted that ae one coefficient 
of velocity distribution enters into the final equation, the one related 
to the distribution of the squares of velocities in a cross section. 

This equation was derived on the basis of the assumption that the 
flow was turbulent, which is the case usually met with in practice. 
Nevertheless, the equation also applies to viscous laminar flow, since 
the dynamic al equations of motion, eq 13, are the same for both types 
of flow, exe ept that in laminar flow the quantities Dery » » +) aFe Fre- 
placed by S.,, ... ,rather than by R,,, .... Its tacitly assumed 
that the flow pattern is such that the assumptions made concerning 
the Reynolds stresses in deriving eq 22, 23, 24 are also true for the 
viscous stresses. Finally, in laminar flow a ‘and \ become quantities 
that may be calculated rather than ones that must be experimentally 
determined. 


2. COEFFICIENT OF RESISTANCE AS A FUNCTION OF H 


If it is assumed that the coefficient 4 defined by eq 45 is a function 
of H only, and is therefore independent of 0U/dt and dU/dz, then \ can 
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be evaluated from the ordinary parameters of uniform flow. Since, 
when the flow is uniform, U and H are constants, we have from eq 47 


. AU? 
— Th,” 
Urs _ [2 

Voi, VX 


where U, and R,, respectively, are the mean velocity and hydraulic 
radius, for uniform flow. 

We assume further that the propagation of the waves is restricted 
to that type of turbulent flow in which the friction is determined by 
the surface asperities independently of the Reynolds number. The 
resistance law is then represented adequately by Manning’s formula 
in the following form [5, eq 61, 80]: 


ne 217)” 


(48) 


where k, is the equivalent sand roughness [5, p. 713] of the channel 
surfaces. 
Equating the right-hand members of eq 48 and 49, we obtain 


» \1/3, 
A= 0.0303( F) 


Differentiating with respect to Hy, the value of H for undisturbed 
flow, we have 
ad» ~~ r dR, r() 
dH, 3R, dH’ © 


since R, is a function of Hy. For a wide rectangular channel, R=H; 


and 
—_— 7 is 
dH, 3H ie 


3. VELOCITY OF PROPAGATION OF SHORT WAVES WITH 
NEGLIGIBLE CURVATURE 


Assuming an initial condition of uniform flow in a very wide 
rectangular channel, let Hy be the depth of the liquid and U, be the 
mean velocity of flow in a cross-section. If the flow is disturbed, this 
can be expressed by writing: 


H(z, t)= oth(a, t), 
and (52) 


Uta, t) — U)+u(z, t), 


for the depth and for the mean velocity. The quantity wu will, in 
general, henceforth signify the variation from U, the undisturbed 
value of U, and h the variation from H,, the undisturbed value of H. 

We replace H and U in eq 47 by their values from eq 52. Since the 
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channel is rectangular, A is a linear function of H; since it is w ide, 


R=H; and Hy and Uy are constants, hence eq 47 and 21 become 
( -b oh e x [ oT 2 
gn 40 r+a(l ot ue > ii~ a) - a? ave (03) 


(H+ 4, on Let . (54) 

We now introduce two restrictions of a different kind. First, we 

assume that @ is very nearly equal to unity and hence the term con- 
taining 1—a@ may be neglected. Secondly, we assume that 


Ou Ou ‘aie 
(I) us <5}? (55) 
. O* Se 
(IT) I — a7 Soe? 
oh poe dh, 
or ' Ox Ot 
On this basis eq 53 and 54 reduce to 


_ OU, 77 OU | 
Tar U Ar 


(ITT) u 


oh 
oF + Hoot + Uae =0 (59) 


We differentiate eq 58 with respect to z and eliminate the terms 
containing wu by means of eq 59. This gives 


oh 
eta U e+ (Up gh) i 0. (60) 


A particular solution of this is 
h=F(x—«pl), (61) 
where J’ is an arbitrary function, provided w, is a root of 
wo? — 2 Ugg + Uo? —gHy=0 
The two roots are 
w= Uy+v 9th, (62) 


and it is evident upon inspecting eq 61 that the wave is propagated 
——— change of form with the celerity, or velocity of propagation, 

in the positive z-direction, since at time t+ (Az/wo), the wave profile 
is is displaced by Az from its position at time t. When an observer is 
moving with the current, that i is, moving with the velocity Uo, a wave 
may be moving either downstream or upstream with respect to _ 
observer. In the first case, Up<wo, and the wave is descending; i 
the second case, Uy>>wo, and the wave is ascending. Ac noadieay, 
the larger root of wo gives the celerity of a iimenine wave and the 
smaller root the celerity of an ascending wave. In all formulas 
containing an ambiguous sign the upper sign will pertain to descending 
waves and the lower sign to ascending waves. 
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It is now a simple matter to determine the function uw. From eq 
61, we have 


oh ol ‘ 
iia cae eee ' 
ot ° Or 
I 
Substituting this value of Oh/dt in eq 59, it becomes 
rus ( i —— Wy oh wut 
Or H, x ; 
whose integral is ) 
y+ Vo~ Mop 4 #74) = 0 
|: ed 
Since for any time t, at h=0, w=0, then since the liquid is undisturbed 
at infinity /(t)=0, and 
(w»— Uo) } (64 
\~OD ( 
“Hy 
or in view of eq 62 
cart ; 
u=++79gHi» 65 
VI "Fy ( 
which states that the excess of the mean velocity of the particles in a 
cross section under the wave is proportional to the wave height h. 
It will be noted that the above equation for the celerity of waves, 
eq 62, is the same as that developed for the motion of irrotational 
waves in still water [1, eq 56], except that a constant velocity of flow, 
. . . . ie 
(9, is superposed on the wave velocities relative to the water, yg/). 
This would appear to indicate that the friction due to turbulence t 
which is inherent in the undisturbed flow, U,, has no effect on the { 
propagation of waves. In general, the effect of turbulence cannot be 
neglected, but in this case the assumptions made in deriving eq 62 
restrict the size of the wave in suck a way that the effect of turbulence t 
Is negligible. ’ 
The physical significance of these mathematical assumptions will 
now be investigated. Consider first 
Oh , ,Ou Oh 
(IIT) u— ih <—- 
Or or “Ot 
Replacing uw by its value from eq 64, and Oh/dt by its value from eq 63’ 
we have 
ho wo " | 
: «<;5( =. eS (66) 
FA, _ 3 0 &o 


Consider secondly 


u ~~ 
or 


Ou 
ot 


Using the same relations as before, we have 


(I) < 


a. 


Hy, ~Us— w 


Obviously, the first condition is sufficient for the second. 
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Both of these restrictions have no bearing on the fact that the 
undisturbed flow is controlled by friction. ¥ hey are necessary but 
are not sufficient conditions for a wave to move without change of 
form. A more significant restriction arises from the assumption “that 
the right-hand member of eq 53 is negligible: 


(U0, +u)? Ou 
S(Hy+h) Sat 


Expressing the left-hand side in terms of u/U» and h/Hp, and replacing 
\ by its value from eq 48, we obtain, since Ry>=Ho, 


h( uH, (uH, 
oi ma} -(7 93 ou 


Sor" 
lta r 


dI gi— 


Eliminating uw via eq 64, and transforming the right-hand side by 
means of eq 63, we have 


gitt(3- (G1) Nee dh. 
— voles) (14 fr ) HJ 


In contrast to conditions (1) and (II), where we obtained an 
inequality comparing h/H») with a constant term which is determined 
by the characteristics of the channel, we find in this case that both 
sides of the inequality are functions of h. We are also confronted with 
the fact that O0h/dzx generally vanishes at the apex of the wave profile, 
and hence the inequality cannot be satisfied at that point even if the 
left-hand member also vanishes. 

In order to arrive at a workable criterion for determining whether 
the theory is applicable, we integrate both sides of the inequality 
with respect to x, thus obtaining an average condition: 


igH, 27 — ; a lanl, be 
— wo (wy— Uo) , ‘Cr dr «< {ip lr. (67) 


It is not feasible to integrate the left-hand side of this inequality unless 
the function h(z,t) is given. It is therefore necessary to assume that 


2a 


7 Wy (2Q 
(IV) '; 1) ye at Mi (68) 


and eq 67 reduces to 


lh | Oh! 
Bf H,, dn = dr 
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where 


2 wy— 3U, 
B= GHWT (a — a Is) 

If the wave profile is such that (IV) and (V) are not both satisfied. 
then eq 70 does not apply, and the integral on the left of eq 67 must 
be evaluated. Assuming that the wave profile has a unique maxi- 
mum heighth», that 0h/or is everywhere finite, that Z is the length 
of the wave, and that A is always positive, then eq 70 becomes 


Oe 

Hy S5B’ 
where 7 is a fraction such that the volume of the wave per unit 0 
of channel above the undisturbed surface is jh L. This last inequalit 
is the desired average condition, one which restricts the length of thy 1e 
wave. Waves fulfilling this condition will be referred to as shori 
WAVES. 

In addition to this restriction on the length of the wave, we still 
have the three conditions on the maximum height of the wave, 
expressed by eq 66, 68 and 69; the first being derived from the ori- 
ginal restriction (III), and the other two from the approximations 
used in obtaining eq 71. Each of the first two restrictions states 
that the maximum height of the wave must be small with respect to 
a constant determined by the channel parameters. For convenience 
we shall express these constants in terms of U)/aH). 

For descending waves, w)=U,)+-¥gH», and eq 66, 68, 69 become 


abhi h 
(IITa) > B spt) 


| U, 
(IV: ee (==- 2), 
ee) glo 


ie h 
(V) H,<- 


7 ° r (TT 
For ascending waves, w)=U,)—~+gHp, 
and the same three equations become 


®. ot 
ab) <TH 1) 


h_ Uf Up 
7 ) 
nsideas He<gHA .V art?) 


(V) = <1. 
0 


These relations are made more intelligible by means of figures 6 
and 7. The quantity h/H, must be small with respect to all of the 
quantities plotted in the figure, that is, small with respect to the 
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least of them. It will be noted that the theory breaks down com- 
pletely for ascending waves in water of critical depth, that is, it does 
not apply to the case of the hydraulic jump. The fact that it breaks 


a 


Ila Va 





Vy 
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FIGURE 6.—Limitations on height of short descending waves moving without change 


of form. 
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Fiaure 7.—Limitations on height of short ascending waves moving without change 


of form. 


down for U>/-/gH)=2 when the wave is descending is due merely to 
the approximation used in deriving eq 71, and is not inherent in the 
theory. 
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Our three original assumptions thus indicate that if a wave is short 
and of small height, and moving in a wide rectangular channel where 
the undisturbed current is controlled by friction, the waves move 
without change of form, and their velocity is constant. It was also 
assumed in deriving eq 47 that the surface curvature is negligible, 
and this restriction is still in force. The simple theory thus gradually 
loses its applicability as the waves become larger and larger. This 
case will be treated in the later sections. 


4. PROPAGATION OF DISCONTINUITIES OF SLOPE 


One of the simplest types of wave profiles which we can imagine 
superposed on the initial uniform flow in a wide rectangular channel 
is one which is polygonal in form, such as a triangle. We shall inves- 
tigate the motion at a vertex of such a polygon, where the slope of 
the wave is discontinuous. 

We assume that at t=f), the discontinuity of the slope Oh/dz is 
located at the point =z, and the water surface is represe nted by 
two straight lines (see fig. 8). The mean velocity, U, and the depth, 
H, will each be described by two distinct functions; downstream from 
Zo, we write U=U, and H=H,, and upstream from Yo, We write 
U=U, and H=H,. We assume that both H, and H; have negligible 
curvature, and that the channel is rectangular. Under these condi- 
tions, eq 47 applies to each pair of functions, and we write, putting 
a equal to unity: 


oH, , ou, _ . AU? 


+t Us “9 FH,’ (72 


ae 


2Us_,,_ Ue 78) 
9 3H, a 


ee 
= ii 4 


Mil pts I 774 i, f 2 
FicurE 8.—Motion of a discontinuity in surface slope. 


In a similar manner, we obtain from eq 21, since A is proportional 
to H, 


oH, , ,,oHh . - 
Ot = +H =0, 


oH 4 y 2H Heh ‘ 
* i 


4t- 
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Following Hadamard [13, p. 83], we introduce the following nota- 
tion for the change in value of a derivative as the discontinuity is 
traversed : 


Or Or 


tor? or Or Or 


OH, dH, = OU, dU; _ = 


Subtracting eq 73 from eq 72 and eq 75 from eq 74, we have, since 
H. U. and X are continuous, 


oH o0U], ,foU (76) 
3, } ot } LB: } " 
oH], ,faH], , fou 
a - F BS | °; 


where H and U are evaluated at r=, t=t. Since four disconti- 
nuities are involved, we need two more relations. 
In figure 8 the displacement of the discontinuity during the time 


interval At is indicated. If AU, and AU, are the changes in the 
mean velocities, then 


aU = art at (78) 
= 


Al a or 


Uy 
Ar+ ot At, 


Similarly, if AH, and ASH, are the changes in the height, then 


dH, , dH, 


AH, = Ar+ at Aft; 


~ Ox 


oH, , oH, 


4H,.= Oz Ar+ dt At- 


Subtracting eq 79 from 78 and 81 from 80, 


ou ou 


oH oH 
o—| 3 +5) At, 


since U and H are continuous. We denote the velocity of propagation 
of the discontinuity, Az/At, by 2. Thus the last two equations become 


oU7 , au _ (82) 
g a or |= 


oH] , [oH] _ (83) 
Ls tLe J-° 
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We now have a system of four homogeneous linear equations, eq 
76, 77, 82, 83, in the four unknown discontinuities. Since the dis- 
continuities do not vanish, the determinant of the coefficients must 
vanish, and this relation gives 


2 — U= + VoH. (84 } [ 


The velocity at the front of a descending wave, as shown in figure 
9, will now be considered. In this case, 2 is positive, and we see from 
eq 83 that [0/7/07] and [OH/dt] have opposite signs, and since in front 
of the wave, 0H,,/ot=0, 0H,/0x=0, then behind the front of the wave, 
0H2/dt>0, 0H,/or< 0. Now, U=U, and H=A, at the discontinuity, 
and the velocity of the wave front is 


9=U,+79H; (85 


the positive sign is selected, since 2>0, regardless of the magnitude 
of Uy and Hy. We thus see from eq 62 that at the front of such a f 
wave Q=a . This analysis fails to give any information concerning 
the change in the magnitude of the discontinuity. Such information 
may be obtained by using Boussinesq’s equation for the velocity of 
propagation of a wave-volume element. <A typical case is discussed 
in section 1V-6, page 498. 


V 


—— 





MANS 
Figure 9.—The front of a descending wave. 


5. EFFECT OF VELOCITY DISTRIBUTION ON THE CELERITY OF 
SHORT WAVES 


III-3, page 478, it was assumed that 1—a is a small fraction that 
may be neglected. This is generally true in practice, but we now 
drop this restriction and thereby consider the effect of the velocity 
distribution on the celerity of chet waves. Retaining a and using 
the assumptions of section III-3, eq 55, 56, 69, the equation of motion, } 
eq 53 now becomes | 


dh du, 7,du Urtudh_ ” 
Isnt t talos, + (l-a)—p— FY mbes (86) 


( 
( 
In the derivation of the celerity of short waves, given in section ( 


Even though 1—a is not negligible, it is still true that in practice 
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is of the same order of magnitude as unity. From eq 69 we 
have therefore 


w—U, dh h ~w—U) Oh 


i“ ont 2 a 


ng the relation between h and u obtained in the previous analysis, 
eq 6d, this inequality reduces to 


uu oh ou 


eae a ae’ 


ind eq 86 thus becomes 
Oh , Ou . Ou Uo 0 Oh _ 
Jax t ag t erg, + on ot 


Differentiating with respect to z, assuming that a does not vary 
with z, and eliminating the terms containing uw by means of eq 59, we 
find 


Oh 
orot 
A particular solution of this is, as before, 
h=F(r—«ayt), 
where F' is an arbitrary function, provided w, is one of the roots of 
we? —2aU pay + (aU,y?—gHp) =0. 


Oh 


ate all, 0 + (aU? 9gHo) a= 


The two roots 
w= (1-9) Ut VgHo+ (n+1°) Ue, (87) 
where n=a—1. (88) 
This is the formula giving the celerity of short waves when the 
effect of the velocity distribution is considered. It will be seen that 
n is the excess over unity of the ratio of the mean of the velocities 
squared, taken over the different points of a section, to the square of 
the mean velocity. A corresponding expression is given by Bous- 
sinesq [3, p. 285, eq 265], which differs from the above only in the 
values of the coefficients. 
It is preferable that the value of n for a given case be ascertained by 
observation. If this is not feasible, a value may be obtained by 
considering the relation which exists between 7 and Uj and u,, the 
latter being the maximum velocity occurring at the su rface of a wide 
channel, A relation of this kind can be derived if it is assumed that 
in a channel of great width the ‘‘velocity defect,” that is, the dif- 
ference between the maximum velocity Um and the velocity u at a point 
whose height above the bottom is z, depends only on the so-called 
“shear velocity” [5, p. 709] and the relative depth ‘of the water, and 
is independent of the character of the asperities on the channel 
bottom. Mathematically speaking, 


Um—U=Uz f (§), (89) 
where ¢=2/H,, and uy is the shear velocity, 
Uy =VigH. (90) 
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f Multiplying both sides of eq 89 by dé, integrating from z=0 to 
z=H, that is, from 0 to 1, putting 


1 
f, S(Hdg=M,, 


and recalling that the mean value of wis Up, we obtain 
Um— Us= Nig. (91 


Squaring both sides of eq 89, we obtain in a similar manner, using 
eq 88, 31, 26, 


Um? —2UmU, ot (1+ 7) )U;,?= Nou’, 
where 


J . . 
J f(g)*dé=N>. Eliminating the shear velocity, u,, between eq 91 | 
Uv 


and 92, and solving for 7, we obtain 


o=(Na— NEY) oe 


which is the desired relation expressing 7 in terms of u,, and U5. 

To determine the numerical values of N, and Nz we make use of 
some of Nikuradse’s experimental! results concerning the velocity dis- 
tribution in sand-coated circular pipes [14, p. 18-19]. The data in 
table 1 were obtained by selecting three determinations of the velocity 
distribution for each of six relative roughnesses. In Nikuradse’s 
paper the data for each relative roughness are given in a separate 
table, and the three determinations on the extreme right of each table 
were selec ted, except in Zahlentafel (table) 13, where the first one and 
the last two were chosen. In each case, velocity-defect ratios were 
computed for each value of y/r, which corresponds to z/H, in our 
notation, and the mean of these 18 determinations was entered in 
table 1. Plotting curves for /(¢) and [f(¢)]* from these data and in- 
tegrating numerically, it is found that in circular pipes, N,=2.52 and 
N,= 12.29. Since the ‘ ‘velocity defect” relation is the same both for 
circular pipes and for open channels of great width these values of 
N, and N, may be substituted in eq 93, which becomes 


2 
n=0.94 ur 7 (94) 


An expression for 7 in terms of \ may also be derived. Eliminating 
Um between eq 91 and 92, we have 


. uU 
- jee (N2.— N,’) ae 
Eliminating u, and U, via eq 90 and 48: 


n= (M—i)5) 
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which becomes, when we replace N; and N, by their numerical values 
n=2.972. (95) 


TaBLE 1.—Velocity defect as a function of distance from channel bottom. 
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It will be recalled that in deriving eq 62 of section I-3 it was 
assumed that a is nearly unity, that is, that » is small, and that this 
restriction was removed in deriving eq 87. By comparing these two 
equations we can establish an approximate relation which will show 
how small » must be in order that eq 62 may be used. . 

Equation 87 may be approximated by 


, ' r =e nn Us 
=(1+4+n)l ot Vg 1+ by gH)” 


and this reduces to eq 62 if 


4+? U,? — 
nl t+ ml Jol <U,+ VgH. 
a is 0 


Dividing by Up, neglecting 7? with respect to 7, using eq 48, this 
becomes 


{1 tay i i +. Very (96) 


or, replacing 7 by its value from eq 95, 


en 1 /2i Ny " 
2.972 Itsy5 “i+ 7 (97) 


Bazin’s inexhaustible store of experimental data furnishes an 
example [15] of the improvement obtained by using eq 87 instead 
of eq 62. The observed values, together with the results of the two 
formulas, eq 62 and eq 87 are shown in table 2. These experiments 
were conducted in the discharge channel of the Grosbois reservoir. 
This channel was constructed of masonry, with the steep slopes of 
0.101 in the upper part and 0.037 in the lower part. Before entering 
the channel, the water drawn from the reservoir flowed over an apron 
11.20 m wide just below the gates, and which narrowed down to 1.80 
m (the width of the channel). The axis of the apron deseribed an 
arc of 90°, with a radius of 40.50 m, in order to join with the channel. 
As a result of the narrowing of the cross section and the change of direc- 
tion of the apron, the flowing water was subjected to various reflec- 
tions at the walls and these complex movements set up waves which 
passed rapidly down the channel at fairly regular intervals of about 2 

524586-—42-—_5 
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sec (31 per minute). These waves were rather small in height, but 
their transit was easily observed owing to the slight noise which 
accompanied them. 

The observed values of w) given in the table 2 are the averages of 
six individual observations, the + term representing the standard 
deviation. The values of 7 are computed from eq 94. The table 
shows that the values of w computed by the simple formula, eq 62, 
are considerably in error, whereas the values of w) obtained from 
eq 87 compare reasonably well with the experimental values. 


TABLE 2.—Erxperimental and calculated values of wo for torrential flow 








Experimental data (g=9.81 m/sec 2) 


Um 


m/sec m/sec 
6.25 +. 07 
4.32 +. 08 
5.75 +. 06 


Error 


y m/sec % 
—22.9 195 6. 63 | +6. 1 
—8.3 07 4. 36 +0.9 
—13.0 : } 7 








IV. SECOND-ORDER THEORY OF WAVE PROPAGATION 


1. EQUATION OF MEAN FLOW WITH APPRECIABLE VERTICAL 
ACCELERATION 


Assuming that the liquid is flowing in a wide rectangular channel 
(v=0, and R,,=90), that the variation with z of the z-component, 1, 
of the local velocity is negligible, or u= U(z, t) and that the assump- 
tions made in section III—1 regarding the Reynolds stresses and the 
orientation of the axes still hold, then the Reynolds equations for 
turbulent flow, eq 10 become 


OU, ,0U_.._ 1 O(p—R,.) , 1 OR, , 
“Ot +l or I> es Se” (98) 


o— OP —F:2) 
Oy . 


and 
ow, 7,0w, dw _1 O(p—R,,) sie 


pea i xsi Nae en 
Ot ° ole p Oz (99) 


We shall eliminate p from the eq 98 with the aid of eq 99, and then 
eliminate z by integrating over the cross section. 
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We begin by evaluating w in terms of the other variables. Since 
u=U and v=0, the condition of continuity eq 8 reduces to 


ow ou 
Oz Or 

Integrating with respect to 2, we have, since w=0 when z=0 
ou 

2—) 
Ox 

which, using a condition of continuity of mean flow, derived from eq 21 

on the basis of the assumption that A is proportional to H, 


OH, ll ou 


cade oH 1 yok 
H\, ot Or 
Replacing w in eq 99 by its value from eq 101 and neglecting squares 
and products of the first derivatives of H and U compared with the 
second derivative, we obtain, since H and U are functions of z and 
t only, 


’ 


w= — 


reduces to 


(101) 


O(p—R,,) _ _ p2U?(OH , 2 OH, 1 OF 
02 Nat TU aratt T oF ) 


Multiplying by dz and integrating, 


ay ee 2 eH, 1 oF 
H? ( 


p—Res=Pa+ eg(H—2) + 0p a tT dott T? oF) 

since at 2=H, p—R,,=p,. We now differentiate the above equation 
with respect to x, again neglecting the products of the lower deriva- 
tives as compared with the higher derivatives, and the value thus 
obtained for 0(p—R,,)/Oz is inserted in eq 98. The resulting equation 
is multiplied by dz and integrated from z=0 to z=H. Since (R,,),-0 
equals 79, the shear on the channel bottom, and (P,,),-4 is zero, we get 


OU *" pou oe, fUsee,.2 On ,1 CH .. 
ot or art 3 Vor tO arott UT orae)—"— 2H’ 
upon replacing 7) by its value from eq 45. 

This is the equation of mean flow for a wide rectangular channel 
where the vertical accelerations are not negligible but the horizontal 
velocities in a section are practically uniform. 


(102) 


2. HEIGHT OF WAVE AS A FUNCTION OF x AND t¢t 


Replacing 0U’/dz in eq 102 by its value obtained by solving eq 100, 
and multiplying through by H, we obtain 


OU ,0H_, oH , PU eH 
Hy U apt GH- Uy, + 3 
23H , 1 3 


H ere 
+ iota t Waren) =H ZU aia 
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We again introduce the expressions 


H(r, t}=Ho+h(a, bt), 


and 7 ; 
U(z, ) =U,+u(z, 2). 


Substituting in eq 100, we have 


Oh , ;, Ou Ooh , ,ou 
Ot | fs, te Ou | hse a 


while eq 103 becomes 


Ou Oh, , sa Oh , , Ou yor. - oh 
0 ot 7 Ot 7 (gH,—l 0 Jay hoy a +- (gh— 2l t)e, 


H 
HU (orh Oh 4 7 orh " 2s 
wd 3 ox +O ae ort soon) -gll— (105 


when w? is neglected with respect to U,?, and in the term containing the 
third derivatives, u and A are neglected in comparison with U and //. 
Various assumptions as to the nature of the right-hand member will 
be made in later sections. 

The terms of small order in the left-hand members of the above 
equations, that is, the last two in eq 104 and the last four in eq 105, 
are simplified by introducing values from the first approximate solu- 


tion: 


=— iy, 
yy, (64 


FH, 


and 
oh oh 


ot Or 
From the latter and eq 61, we obtain the further relation 
oh 2h 
ot? ie) Ox? 
For the sake of brevity, we introduce the new variables 
3U, 


Oo 


(63) 


ky =|]— 
and ; 
"9 —— 1 -_ l * 


Wo 
and we also make use of the equation 
gHy= (w— U,)?, (107) 


which is derived from eq 62. On the basis of the above relations, then 
eq 104 and 105 become 
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oh ou oh 
57 tise t Us” + (wy— oe (a )= 0: (108) 


: \o k, h? 
Het — Uer+ Hy—T aye = + 6p (629 — o) 2 3 Th 


—— 20) ‘Tp 


(109) 
We now differentiate eq 108 with respect to ¢t. In differentiating 
the last term, which is of the second order, we use the relation 


ry: ( > 

) 110 

(3)= ~~ ss (110) 

since these small terms may be treated as functions of r-—wot [1, p. 73]. 

We next differentiate eq 109 with respect to z, and subtract this equa- 
tion from the one obtained from eq 108. The result is 


O° _ Oh rok 
pt? Doan ot — GMy— U0") 52 


kH,? 07h re 
~ngihaigon EIS IZA 7 “(a +" : 3 )=— 5( ti nU®)s (111) 


which is the fundamental equation for determining / as a function of 
r and ¢. 

This equation, which was derived from eq 102, has the same range 
of validity, namely, it applies to flow in wide rectangular channels 
where the horizontal velocities in a section are practically uniform. 
This equation takes account of vertical accelerations, finite wave 
heights, surface curvature, and friction. 


3. CELERITY OF A WAVE-VOLUME ELEMENT 


In order to integrate eq 111, the concept of celerity of a wave-volume 
element, an essential simplification due to Boussinesq [3, p. 451], is 
introduced. To define this concept, we consider the volume of the 
wave per unit width of channel, above the undisturbed primitive level 
of the liquid and between z, an arbitrary point, and ~, for a descending 
wave, and between z and —o, for an ascending wave. Denoting 
this volume by o, we have 


o= [hie and =|. hdz, 


for descending and ascending waves, respectively. Recalling the 
convention of signs introduced in section JII-3, the above equations 
become 


+o 
eu +f hdz. (112) 
Zz 


Imagine a plane, normal to the channel axis at the point 7, and which 
moves in such a way that the wave volume in front of it, 7, remains 
constant. This velocity is the celerity of the volume element oc, and 
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is denoted by w=wo(z,t). Differentiating eq 112 with respect to ¢, 


we have ‘ 
+02) 
o=+ f Ode he, 


since by definition dz/dt=w, o remains constant, and we assume that 
the surface of the liquid is undisturbed at infinity, that is, h vanishes 
at +o. Differentiating the above equation with respect to z, we 
have, regardless of choice of es 


oh 
a2 
ot ' or 
which gives the relation obtaining between w and h. 
We now introduce the value of 0A/dt from eq 113 into the first two 
terms of eq 111, thus obtaining 


O° _ 3? aed Oh 
~sr01™ —2l 09,2 dw) — (gH,— Up hod 


| _ BHe ot a 
— ele Wd 5 Ht 3 on) ad GHE-BU") 


This equation is immediately integrable once with respect to x. The 
arbitrary function of integration vanishes as a result of the previous 
assumption that the surface of the liquid is undisturbed at infinity. 
We thus have 


is: , oO 
5p ho) +2 U5. (hw) + (gH)— 


(hw) =0, (113) 


122k 
Ut) >, 


0/2+hy WM, keHT,? Vh\_ ; 
— Wp (wo— Us) 0) s(75 mej 3 er oa )= gl —3U? (114) 


We next transform eq 114 into a form which can be integrated. 
In this process it is necessary to operate only on the first three terms 
in the left-hand member. If we assume that w/w is near to unity, 
we may write as an approximation 


w= w(1+6(z,t)). (115) 
Substituting this value of w in the terms under consideration, we obtain 


cg (h-+ hd) +2 1 2 wy U, Sr Duets hé) + (gHo— U; 22%. 


Now using eq 113 and 107 we obtain 


o o «o ; oh 
— os: (hw) + o>} (hb) + 2a9L Sr (A+hd) + (co? — 229 U5) = 


Applying eq 115 again to the first term and eq 110 to the second term 
(which can be done because hé is a small quantity), and cancelling 
terms, we obtain 


) -o 
ant 2" (hd) + 2 Sz (hd) ’ 





Translation Waves in Open Channels 


which can be written 


5. 
0) 5, (248): 


= Wo (w = l 


Hence eq 114 transforms into 


a) 
I koF i 

Q+kh fh »H,? O7h So 
("y,)+45> Se) sie" 


— FJ 
Wy \ Wy l ‘y) 57 (2M) + Wp (wy — ‘a 


which, when integrated, gives 


24 ky h k.H? oh 1 , . od ») 
oe T+ 6h So ET, (oH BU dx, 


and hence from eq 115 


Oth hk, kyHy? dh a 
a ae: Sat heed, (gtti—3U*)de. (116) 








This is a general expression for determining the effects of wave 
height, curvature of the wave profile, and friction, on the celerity of a 
wave-volume element. The integral containing the friction term 
requires further consideration before the actual wave profile can be 
determined as a function of time. 

It remains to evaluate that part u of the velocity U under the wave 
ineq 52. Transforming the equation of continuity, eq 100, by means 
of eq 52, 

OF +2 (Uo+u) (Hoth) =0. 


The result of subtracting eq 113 from the above is 


out o+u) (Hy +h) —ho]=0. 


Integrating, and using the boundary condition at r= @, for all times 
t, h=0, and u=0, we have 


(Uptu) (Hy +h) —hw= UH, 


ecm Jy 


4. DEFORMATION OF THE WAVE PROFILE, ASSUMING 
gHi-*U%=0 


A natural assumption, in proceeding from eq 116, is that the right- 
hand member (gHi—\U?/2) is negligible. From eq 56 this implies that 
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the waves are short, as described in section III-3. On this basis, eq 
116 becomes 
2+kh, h k.H,? o7h 
wii wf 1 + _ - Ht 6h =) (118) 


or when we introduce the values of ki, 2, and wy from eq 106 and 62, 


: — 3h , H,? ovh 
w—l o— + Vg 1 +aH+ Gh 53) (1 19) 
Comparison with an earlier paper [1, eq 112] shows that this expres- 
sion is equivalent to that for the propagation of an irrotational wave 
relative to the mean velocity of flow, Up. The assumption that 
gHi—rU?/2 is equal to zero is thus tantamount to assuming that the 
only effect of friction and slope is to superpose a constant velocity U, 
on the wave, which itself is propagated without any change of shape 
due to frictional damping or to the z-component of gravity. 

The deformation of the wave profile, in a short time interval, is 
determined by using eq 119. Equation 113 can be written 


oh, oh, O 

“Ot ! wos rd (w— wo) =0. 

For small values of ¢, the integral of the above equation is 
h=f(e—cat)—ta{h(w—en)] 120) 
v=} (x— wot) Dl (e— oo) J. (120) 


In verifying this solution, we must use eq 110 when differentiating the 
second term on the right, which necessitates assuming that this term 
is small compared with the first term. In other words, the magni- 
tude of the time interval in which this approximation is permissible 
varies inversely with O[h(w—w,)]/or. Clearly, since at t=0, we have 
h=f(r), the function f represents the wave profile at the fiducial 
instant. Finally, using eq 119 and 62, eq 120 can be thrown into the 
form 


(121) 


| _ aah? , Heo 
h=f(xr— wot) F ty gs srr oe s) 


For a discussion of the deformation of a wave of negligible curvature, 
propagated without damping, the reader is referred to the earlier paper 
[1, p. 89]. All the formulas on the page cited are directly applicable, 
provided the expression ¥gH(1-++ ...) isreplaced by Ujt+VgM(1+...), 
wherever it occurs. This, again, is merely equivalent to superpos- 
ing the velocity of flow Up, on the velocity of propagation of the 
wave-volume element. 


5. EFFECT OF FRICTION AND SLOPE OF BED ON THE PROPAGATION 
OF WAVES 


In the above solution of eq 116, the term gHi—U?/2 was ignored. 
If the wave is short, the influence of this term is negligible, as explained 
in section III-3. For long waves, however, these effects become 
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appreciable at distances sufficiently removed from the head of the 
wave. The problem of determining these effects is considered sepa- 
rately for the two cases, (I) 7>0 and U)>0, and (II) i=0 and U,=0. 
The first case will now be considered. 

From eq 52, we have, after neglecting u? in comparison with U,?, 


gHi—21 gH — S47 oo tghi— dU ou. (122) 


Considering the first two terms of the Taylor expansion of \= d(H) 


at Ab, aN 
A= Not (a) (123) 


where \=A(Z{), and dd/dH is evaluated at W=H). Replacing X by 


the above value, and neglecting the small hu term, eq 122 becomes 


gln— 7 ?—gHyi- doy ve +gh i} pa Uy? — ry Upu. 


~ 


From eq 48 we have 
a 
gHi—szl oO, (124) 


2 


that is, in uniform flow the effects of slope and friction are in equilib- 
rium, and thus 


dx 


gHi—2U®=ghi—J pyphUg— NU 


Evaluating d\/dH by means of eq 51 and 124, and introducing the 
value of u from the first approximation, eq 64, the above equation 

becomes 
eo “ae es ™ 
gm—sU'= —2(wo— Up) hpi, (125) 

where 
3a,—5U, 

—<- —_ ri Qt. (126) 
ji 3l g(Wo— L 5) 


Since from eq 112, do/dr= +h, eq 125 can be written 


- Ary me Fe ft 
gHi-2U?= £20 (a — Us) (“), 


Substituting this in eq 116, we have 


o=a(1 — +h, h ko)? Oh hy ) 


4 Hy" 6h Ox?" woh 


or introducing the values of k; and k, from eq 106, and using 


—_ _ 3h He eh\ _ oft). 
w— U,= +(voH(1 + oH, t Gh 51) — 3 ) 
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This integral of eq 116 is due to Boussinesgq [3, p. 451], except that the 
value of fs given above is new, being based on Manning’s quadratic 
law of rE AT tae developed in section ITI-2. 

We now consider the sign of the quantity of,/A. If the intumescence 
is everywhere positive or everywhere negative; that is, for all z, h=>0 
or h<0, then o has the same sign as h, and o/h is positiv e. Atthe head 
of the wave this ratio vanishes, since ¢ is of higher order than A, and it 
increases toward the rear of the wave. Since f, is constant for a given 
channel and depth of water in uniform flow, the absolute value of the 
term varies in the manner just described, and its sign is the same as that 
of f;. Upon examining eq 126 we find that f, is positive if (1) Uj >w 
or (2) Ujp<wo and 3w)>5U5, and negative if (3) Ujp<wp and 3a)<5U’, 
Comparison with eq 62 shows that (1) applies to all descending waves, 
(2) to descending waves where U,’< (9/4)gHp, (3) to descending waves 
where U,?>(9/4)gH,. Thus, except in the last case, the combined 
effect of friction and slope on the damping of the wave is to cause a 
diminution of the absolute value of the velocity of propagation with 
respect to the velocity of flow of the various parts of the intumescence, 
the effect increasing with increasing distance from the head of the 
wave. The last case is of special interest. It has been shown else- 
where [16] that if U,?>(9/4)gHo, the waves become unstable, a 
phenomenon which is apparently closely related to the origin of roll 
waves. 


6. DEFORMATION OF A STRAIGHT SLOPING FRONT 


We suppose that the wave is descending (fig. 10). Denoting the 
slope of the front by 8, B=—tan ¢g. Since we are only considering 
waves whose profiles lie above the free surface, B is always positive. 
Denoting, at t=0, the position of the front of the wave by 2, and the 
value of the slope by Bo, the equation of the wave profile is thus, 


h=B)(1%—2), LEX, i=0; 
h=0, 12%, t=0. 


Equation 112 then becomes 


=. 


o=0 LZ2=%, f=. 


(130) 


Introducing the equation for the velocity of propagation of a wave 
volume element, eq 128, we obtain from the above equations 


oma 1+(7 7H 3 )eo—2) } tS%, t=0. 


If we write 
2+k, ti 


A= GH, “Po 3’ 


then the last equation becomes 


W— Wy = 99 (%p— 2), LEN, §=9. 
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FicurE 10.—Deformation of the linear front of a descending wave. 


Now let the position of the front of the wave at {=At be 2)’ (fig. 10). 
The section which at time ¢=0 had the volume o in front of it has 
moved from z to x’, and the height has changed from h to h’. Since 
the wave front moves with the velocity w. (section III-4), while the 
section at x moves with velocity w, we have 


/ 


Io —x’ = (wy—w)At+(%—zr), r1<ay, t=0. 
Using eq 132, we obtain 
Io’ —x’ = (1—H At) (%o—2), t< Xp, t=. (133) 


But the volume ¢ in front of z is equal to the volume « in front of z’, 
and thus we have from eq 130 and 133 


oO} = Po(z,! —z')*(] + 26)At). 


Comparison with eq 130 shows that during the time At the linear 
character of the wave front is preserved, but the slope has changed 
from By to 8’, 


B’ = (1+ 26)At) Bo. (134) 
Putting 


19 1 4B 
B’= Bo+ Fat, 


we have from eq 134 


dp, on 
+, — <0 ==()), ioe 
dt 9080, oy 


which enables us to determine the variation of the slope. Since fp 
is positive, we find that if @,.>0, the wave is becoming steeper, if 
6=0, the slope remains constant, if <0, the wave is flattening out. 
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It is desirable to translate these conditions into relations iny olving 
Bo. 7, and »X. Using eq 106, 126, 62, and 48, eq 131 becomes 


3/9 - 
— 4\ I, 1 (g,— -, (136) 


.: = oe - 
where J=i Gy =), Since g, Ho, and Bp are positive, we thus find 
‘ ray} t 


that if By) >TJ, %>0; if B=T, %=0; and if &<J, 6) <0. Accordingly, 
we find that, 


ef] , A. ; ; 

Bo> (3 d/21— 5) the wave is becoming steeper; 
7. Ok 

Bo= (37 d/2a— = )is the slope remains constant; 


Br<(24/ r/2i-2)i, the wave is flattening out. 


We now derive an expression for 6 as function of time. Since it 
has been demonstrated that the wave front conserves its linear char- 
acter, eq 135 holds for any value of ¢, and we obtain, using eq 135 
and 136 and discarding the subscript which refers to t=0, 


dp 3 
dt Ae — 


The solution of this differential equation for the initial condition 
B=, at t=0 is 





are 


0 


Bol 
Bo— (Bo—L) exp a, It 


If J=0, the solution of eq 137 becomes 


2 ii 


and the slope continually increases. 

We next suppose that the wave is ascending (fig. 11). Denoting 
the slope of the front by 6, B=tan y. We again suppose that the 
wave profile lies above the free surface, which makes 8 positive. 
Using the same method of reasoning as in the above, we obtain cor- 
responding to eq 135 and 136, 


dp_ 
dt —s =26oBo, 


3 
A= ay i —J +f), 


J=Gy +5): 
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Thus, if Bo > =~ 5) the wave is becoming steepe1 


| ey. ys . 
yg Je the siope remains constant; 


/ 


2 ih, &X. , : 
3V arte the wave is flattening out. 





Figure 11.—Deformation of the linear front of an ascending wave. 


The expression for 8 as function of time, becomes 
(140) 


and if J=0, which for ascending waves can only be the case if the 
slope, 7, vanishes, the latter reduces again to 


— a 


3/9, 
“ova 


and the slope of the wave is continually increasing. 
If it is desired to express 6 as a function of the distance 0’ —2 tra- 
versed by the wave, this is obtained by introducing (z9’—2,) W in 
place of ¢ in the various expressions giving 6 as a function of t, since 
the head of wave is traveling with the velocity a. 


7. CHANGE OF HEIGHT OF AN ABRUPT WAVE FRONT 


In a wave having an abrupt wave front there is a discontinuity of 
height at the wave “front (fig. 12). If for some considerable distance 
behind the front, the height is constant, the change of height with 
time, or preferably with the displaceme nt of the front, can be readily 
determined, using the Boussinesq formula for the velocity of prop- 
agation of a volume element. 
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FicurE 12.—Deformation of an abrupt wave front 


It is thus assumed that at t=0, in the initial portion of the wave, 
h is constant, and the wave is descending, that is, 


h =h,, = Lo, t=O; 
A=0, Z2%, t=0. 


Since the wave is descending, eq 112 becomes 
o=ho(m—zr), LX, t=0; 
o=0, 72%, t==0. 


The velocity of propagation of a wave-volume element is thus from 
eq 128, 106, 62 


34 gH, ho , (x—4X)f,; eee deta d . 
w of 1 a ia 8 H,t = >» £S%, t=. (142 


From this it follows that the velocity of the wave front, w,, is 


— BV gH, ho aa (143 
w= 1 = 4 oa EE ’ i=). 


Since the height is constant behind the wave front, 0h/Ox=0, and we 
have 
dh__ oh 
dt ot’ 


144) 
r<a, t=0, ane 


and from eq 113 


Oh _ 
ot 


(145) 


hoo Z< x, $==0, 


Differentiating eq 142 with respect to z, Ow/0a=/f;, and eq 145 now 
becomes 


dh (146 


di J thos S< le cap. 
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We consider two distinct values of x at t=0, both behind the wave 
front, denoting them by 2, and 2, and write hi=h(2;), hg=h(22). 
Af er an interval of time At, let the height of the corresponding wave 
elements be A,’ and h,’. Now 


dh 


b} —h,= dt 


At= —fih At, 


and 


h,’— he =Gat= —fih At. 


Since hi=h,=he, we see that h,’=h.’. This shows that the con- 
stancy of height of the wave front is preserved. Replacing ho by h in 
eq 143, we see that the velocity of the wave front w, is 


oof 1453 ae = b =t, (147) 


Again replacing hy by h, eq 146 now becomes 


a= —f,h, r<avt { w,dt, (148) 


Integrating this equation, and using the initial condition h=ho at 
t=0, we obtain h/ as a function of ¢ 


t 
h=heM, r<ayt | edt (149) 
0 


If it is desired to determine A as a function of the distance 2’ —29 
traveled by the wave front, we write from eq 147 and 148 


dh _—f,h 3 v9 h - 
dry’ 1 4 Wy Ay iat 


since dao’ /dt=w,; and 3hVgH,/4wFhy is required to be small with 


respect to unity. The solution of this differential equation for the 
initial condition h=ho, at 2’ =p Is 


13 Vg h 
h 4 w@& Th , xp [A= (%o— - I’) )] as 
- _—— , (151) 
i oo 3 3 V¥gH f° 


rl Wo 


8. SOLUTION OF THE WAVE EQUATION WHEN SLOPE AND FRICTION 
ARE NOT NEGLIGIBLE 


The significance of introducing w, the celerity of a wave-volume 
element, lies in the fact that the integration of the wave equation, eq 
111 becomes equivalent to meaty: 


Oh 


—? 2 (hws) = (113) 
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where w is given by eq 116, or in this particular case by eq 128. 
Although the general solution is complicated, if it be assumed that 
O*h/ Ox is negligible, compared to 6h/H,’k2, the problem is consider- 
ably simplified. This assumption requires that the curvature of the 
wave profile be small. We thus obtain from eq 129 and eq 62: 


—_ ; 3h 
wh=hUp+ (w—l 1+ 577, hE oh (152) 


and eq 113 becomes 


e)/ : ¥é 3h oh 
Bit | Uot (eo —Ui)( 14377.) Be=—Mi, (153) 


since 0¢/O0r=+Fh. The general solution of this differential equation is 


-, dh e4,—1 
r= («04 (wy— LC DT, —— t-+(he,), (154) 


where ® is an arbitrary function. This integral is due to Boussinesgq 
[3, p. 457, eq 392] and Favre speaks of it as one of the Boussinesq 
integrals [6, p. 106]. 

In determining the arbitrary function ®, two cases are of particular 
interest; (1) the channel is infinite in both directions, and (2) the 
channel has an entrance and is infinite in only one direction. In the 
first case, if for t=0, h is a known function of z, —~<r<o, then & 
can be obtained directly. Let the function be h=W(r), t=0. Put- 
ting t=0 in eq 154, we have x=9®(h), t=0, hence in this case ® is 
simply the inverse of ¥. In the second case, ® is determined from 
the boundary condition h=y(t) at the channel entrance z=0, but in 
this case the evaluation is indirect and tedious. At x=0, eq 154 


gives 
-, oh ex Rent 
(heh) =—| oo + (Ud ggg PY ‘|: 


In general, the form of the equation does not permit solving explicitly 
for either ¢ or A, hence it is necessary to evaluate @ numerically. By 
considering values of ¢ in the range — ~ to + ~, and the corresponding 
values of A at the entrance, using h=y/(t), the magnitude of the 
independent variable can be computed. Corresponding values of ® 
are obtained from the right-hand member of eq 155 by inserting the 
same values of h and t used to evaluate he‘. 


9. INTUMESCENCES PRODUCED BY OPERATING LOCKS IN A CANAL 
CARRYING WATER 


We now consider a practical example of case (2) above which can be 
treated by a simpler method. In many cases the canals leading to 
and from hydroelectric plants are also used for navigation, and a 
lock is provided to bridge the difference in head utilized by the plant. 
Whenever this lock is emptied, a positive wave is produced in the 
downstream canal, and, whenever it is filled, a negative wave is pro- 
duced in the upstream canal. The former, of course, travels down- 
stream, whereas the latter travels upstream. 
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Because of such design considerations as the possibility of over- 
topping the embankments, it is desirable to predict by computation 
the magnitude and deformation of these waves. Ordinarily the 
discharge from the lock is a continuous function of time, the exact 
relation being obtained from computation or model experiments. The 
waves produced by operating the lock are small in height and long 
enough so that the curvature is negligible, and thus eq 154 is applicable. 

As far as this analysis is concerned, the lock constitutes a termina- 
tion of the canal, and the boundary conditions of case (2) apply. 
The complications of this case are eliminated by a method due to 
Favre [6, p. 112], which is based on the fact that the distance traveled 
by the wave is very large compared to its length. The initial shape of 
the wave is computed for the initial time t=0, in the form h=¥(z), 
assuming that the effect of friction is negligible during the short time 
interval required for the formation of the wave. This function is then 
used to determine ®. 

Let the discharge at the lock per unit width of channel be represented 


by the function 
q=f(t), z=0. (156) 
As a result of the above, we may also write 
t=F(q), z=0. (157 
It follows from eq 154 that the form of the wave at the moment 


it is produced and shortly thereafter, neglecting the effect of friction 
during the genesis of the wave, is represented by 


r—| cyt (a Us t+a(h), (158) 


since Lim exp(f,t—1 
ti>0 tf, 


in @ cross section beneath the wave is 


‘1. From eq 52 and 117 the mean velocity 


U= te ihe (159) 


Neglecting the effects of friction and curvature in eq 129, this becomes, 
using eq 62 


U=U)+ (w— (1+ Fr rok (160) 


The discharge per unit width at z=0 is thus 


q=HU=U,(Hy +h) + (o— M(1+a7 PSO, 2=0. 


Comparing this equation with eq 157 we have 


t=F(Us(Ho+h) + (o— (1+ ze YO 2=0, 
524586—43——_6 
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which expresses h as a function of t for z=0, and thus serves to deter- 
mine the arbitrary function @ in eq 158. Placing z=0 in that equa- 
tion, we have 


$(h) =—| «+ (w»— Up) om | (162) 


and from eq 161 


> 7 . 3h ' 
(0) =—[ ont (oo Uidsty [Pollo +) + (oo — Ua 145% Wo. (103) 


Since the function # in eq 158 is the same as the function © in eq 154, 
we write for the second term in the right-hand member of this equation, 


& (he**) )=—| e+ (oo 0 Se TT, \FWolHo + he’) 


+ (wy— Uo\(1 ae Tie (164) 


The solution of a typical problem will now be described in detail. 
It is desired to find the wave profile in the canal at any instant ¢ when 
the discharge is given as a function of time, Q/b=q=f(0), 6 being the 
width of the free surface of the canal. We shall assume that the 
discharge first increases and then decreases, as shown in figure 13. 
This will create a positive surge in the discharge canal and a negative 
surge in the supply canal. The origin in figure 13 is conveniently 
chosen to correspond to the instant of maximum discharge. The 
abscissa ¢t is, of course, expressed in seconds, and the ordinate, q, in 
any convenient unit of discharge per unit width, such as square feet per 
second. In eq 154 and 164, upon which the solution depends, U and 
wy will have to be expressed in corresponding units of velocity, such 
as feet per second. 
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Fic. 13.—Discharge from supply canal as a function of time: Q/b=q=f(t), eq 166 
or t=F(q), eq 157. 


The positive surge in the discharge canal will now be considered. 
The quantities Us, Ho, R, and 7 are determined from measurements of 
the steady flow. The coefficient of friction, Xo, is then determined from 
eq 48. All the quantities needed to determine wy from eq 62, and f, 
from eq 126 are thus available. Note that the positive sign i in eq 62 
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is selected, since the wave will be moving downstream. A small posi- 
tive value of A is now selected, which we shall denote as h;. The argu- 
ment of the function F in eq 164 is computed, using the value of ¢ 
corresponding to the instant for which the wave profile is desired. As 
is seen from eq 157, this gives a corresponding value of g, which we shall 
call q. If q is greater than the maximum value of q which is actually 
obtained, this signifies that A, was higher than any point of the intu- 
mescence. It will be advisable to select h; very small at first and then 
increase it until the top of the curve in figure 13 is attained. The values 
of t corresponding to h, will be denoted by t, and t,, since, there will 
generally be two of these, the first positive and the second negative. 
This is clearly indicated in figure 13. These two values, t;; and t,., are 
thus the desired values of the function F occurring in eq 164, and the 
two values of the double-valued function ® can now be obtained by 
replacing A in the coefficient of F by the quantity h,;. Turning now to 
eq 154, replacing h by h; and t by the value corresponding to the instant 
for which the wave profile is desired, we obtain two values for z, cor- 
responding to the original value h,, which we shall denote by 2 and 
ty. These are the abscissas of the points on the profile having the 
height h, (see fig. 14). The whole process is repeated for greater 
heights, he, hz, . . . , until sufficient points are determined to estab- 
lish the contour of the wave. 

The process for determining the negative surge in the supply chan- 
nel is exactly the same, except that in eq 62 the minus sign is chosen, 
since the wave is traveling upstream. The quantities h,, ... .and 
Zi, T2, - . - , Will be negative, as shown in figure 15. 
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ee 
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FicurE 14.—-Profile of descending wave in discharge canal. 


10. EFFECT OF FRICTION ON WAVES WHEN THE SLOPE OF THE 
BED VANISHES 


We now consider case (II) of section IV—5, 1=0 and U;>=0. The 
determination of the positive sense of the z-axis on the basis of the con- 
ventions in section I[I—1 fails, since there is no flow. Furthermore, 
the terms “ascending” and “descending” waves explained in section 
III-3 fail to have any significance, since we have no method of dis- 
tinguishing the upstream from the downstream direction. We thus 
are forced to make an arbitrary choice of the positive z-direction, and 
consider downstream to be in this direction. A descending wave thus 





508 Journal of Research of the National Bureau of Standards 








. >. /ff) == i 
z=|,}] => 


i= 
MIS (ES\\\\ } 


2N=7aNs — O x 


FicgureE 15.—Profile of ascending wave in supply canal. 


becomes one for which w>0, and an ascending wave one for which 
w)<0, that is to say, a descending wave is propagated in the positive 
z-direction, and an ascending wave is propagated in the opposite 
direction, since our hypothetical observer is stationary and the down- 
stream direction coincides with the arbitrarily chosen positive sense 
of the z-axis. 

Since i=0 and U,=0, it follows immediately from eq 52 that 


gHi—2U2= — Put. (165) 


It has previously been tacitly assumed that U was always positive, 
that is, uw was smaller than U>, so that U and U, had the same sign. 
Since in this case U,=0, this assumption is no longer valid. The 
sense of the friction is always opposite to that of the velocity, and thus 
if wu is negative, the sign of the right-hand member of the above 
equation must be changed, since we wish to keep J positive, regardless 
of the sign of 7». From the first approximation, eq 64, ua h/Ho, and 
thus the velocity must have the same sign as wh. We introduce the 
new parameter 


f=, if wh >0; 
le 
fens if wh<0; 
and obtain from eq 165 
ss (G)s (167) 
g 9 ai 8 —s Hh 2° ] 


The difference between \,u? and dw? is a quantity of the order being 
neglected in this analysis. If we introduce the quantity 


+o 
o’ = + {iids, (168) 





Translation Waves in Open Channels 
then since do’/dxr= +h’, we can write eq 167 in the form 
se sue € 
M—sU*= +2 fe =) (169 
a SE Ae — 


Substituting this in eq 116, we have for the velocity of propagation 
of a volume element 


— (2k fy beHe’ Oh oo _ o'fe 
orl “4 Hy" 6h d2** w—UphHe) 


Since in this case U)=0, we see from eq 106 that k,=k,=1, and thus 


me i on , tg Oh. o' fy See 
wma ta7 toh Or hi) — 


Consider the sign of woo’fo/hH,’?, which equals w9’o’f2/whH,?. Since 
from eq 168, o’ is always positive, it follows that the sign of the above 
is the same as that of f2/w. But from eq 166 f. has the same sign as 
wot, hence the original term is always positive. In accordance with 
the convention of choice of signs for descending and ascending waves, 
the choice in eq 171 is always opposite to that of the sign of w. There- 
fore, the effect of friction on the celerity # when i=0, U,)=0, is to 
diminish its absolute value by an amount which increases with the 
distance from the head of the wave. 


(170) 


11. SOLUTION OF THE WAVE EQUATION WHEN THE SLOPE 
VANISHES BUT FRICTION IS CONSIDERED 


As in section IV-3, the problem is to solve 


oh, O 
+= (hw) =0 
tt Ox t) 0, (113) 
where w is given by eq 116 or, in this particular case, by eq 171. For 
the reasons mentioned in the same section, we assume that 0*h/Oz’ is 
negligible. Equations 171 and 113 thus become 
oa 4p (179 
w= 1+ F ip ’ (172) 
4H, hH? 


and 


Oh , 3h \ Oh __ h\? (173) 
St (1+3, )oe= Sona) 


respectively, since do’/dr=+h?. Assuming that the wave is every- 
where positive or everywhere negative, then f, is constant, and the 
general integral of the above differential equation is 


Frwy 
7 sh log( 1— I il) . 
ama 14 0E Fa +8] 1 


—Faht 
1,” 
where @ is an arbitrary function. This is Boussinesq’s second integral 
(3, p. 454 eq 393]. 
Regarding the determination of this arbitrary function, we consider 
the same two cases as in section 1V-8: (1) the channel is infinite in 
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both directions, and (2) the channel has an entrance which is infinite 
in only one direction. In the first case, if for t=0, h is a known func- 
tion of r, — ><.r< &, then # can be obtained direc tly. Let the fune- 
tion be h=W(r), t=0. Putting t=0 in eq 174 we have r=9(h), t= 
hence in this case ® is simply the inverse of ¥. In the second case, 
is determined from the boundary condition h=y/(t) at the channel 
entrance z=0; but in this case the evaluation is indirect and tedious. 
At r=0, eq 174 gives 


= a log (1 Ti) | 


? i ag I-53 
1— Srey aT —Se ae] 


H,? He 





In general, the form of the equation does not permit solving explicitly 
for either ¢ or h, hence it is ne cessary to ira ® numerically. By 
considering values of ¢ in the range —~ to +, and the corresponding 
value of A at the entrance, using h= y(t), tie magnitude of the inde- 
pendent variable can be computed. Corresponding values of ® are 
obtained from the right-hand member of eq 175 by inserting the same 
values of h and t used to evaluate 


h 


fo Wy 
l oy 


These methods are precisely analogous to those discussed in connection 
% ce oe DPD 
with eq 155. 


12. INTUMESCENCE PRODUCED BY OPERATING LOCKS IN AN 
ORDINARY CANAL 


In contrast to canalized rivers and hydroelectric canals, such as 
discussed in section IV-9, ordinary canals contain still water except 
when disturbances are set up by the operation of the locks. The 
general considerations at the beginning of section IV-9 (through eq 
157) also apply to this case, except that i=0 and U,=0, and hence 
we must consider case (2) of eq 174. 

The discharge per unit width at the lock is again represented by 


q=f@), z=0, (156) 
or by the inverse function 

t=F(q), r=0. (157) 
Neglecting the effect of friction during the formation of the wave 


and shortly thereafter, we obtain from eq 174 


r= anf 1+ oer t+ P(h), (176) 


lo ie t) 


“ in “0 ~ ae 





Translation Waves in Open Channels 511 


For the mean velocity in a cross section beneath the wave, eq 52 


117, 171 give 
ie =o( 1+ an Cee (177) 


since Uj)=0, and friction and curvature are neglected. The discharge 
per unit width at r=0 is, therefore, 


; 3} 
q=HLU =wol(1+ 77 ) _— 


Substituting this value of g in eq 157, we obtain 


3h y 
t=F(w n(1 +7 ) z=0. 
a a ) me al 
&(h)- an op, eo *=9, (180) 
or from eq 179 


on 3h -) 3h 
&(h)=— + OH. F(«oh(1 +am)) (181) 
Since this function ® “ eq 176 is the same as the function ® in eq 


174, we have for the right-hand member of this equation, 


h 


@p Fy = — (4 foo 
1 ht *oH(1— yt) 


But from eq 176 


Seay) “an “Fe) ; ae 


H,? H,? 


Equations 174 and 182 are used to determine the wave profile at 
any instant ¢, as is described in section [V-9. 


13. LIMITATIONS OF BOUSSINESQ SOLUTIONS 


It is necessary to emphasize the fact that the second-order theory, 
which is based on the concept of the velocity of propagation, , of a 
volume element, is subject to certain important restrictions. The 
most general expression for w is given in eq 116. In deriving this 

equation, the important assumption was made that w/w, is nearly 
equal to unity. 

We shall now turn our attention to the specialized equations derived 
from eq 116. These are eq 118, 128, 171. Examination shows that 
In each case, w/wo=1 + (correction terms). The restriction on the 
theory thus requires that these terms shall be small —— to 
unity. The first term, which is common to all equations, is (2+4,)h 
4H,. This merely requires that h<4H,/(2+k;). The second term, 
which expresses the effect of curvature, requires that (07h/dz’), h 
<6/H,k2. We see from this, that (07h/dx?) must vanish to a higher 
order than h. 

The third term, which appears only in eq 128 and 171, expresses the 
effect of friction. These two equations are applicable only when 











512 Journal of Research of the National Bureau of Standards 


of;/woh or o’f2/hH,?, respectively, is small. But o and o’ can be made 
as large as we please merely by considering longer and longer intumes- 
cences, and if the intumescence has a finite volume, h approaches zero 
at the tail whereas o and o’ remain finite. This shows that the term 
under consideration may be large, or even infinite. In view of this 
fact, it may be expected that these solutions will give correct values 
only in the frontal portion of the wave where o/h and o’/h are small, 
It should be noted that this remark applies with equal force to the 
Boussinesq integrals, eq 154 and 174. 
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MATHEMATICAL TABLES 


Attention is invited to a series of publications prepared by the Project for the 
Computation of Mathematical Tables conducted by the Federal Works Agency, 
Work Projects Administration for the City of New York under the sponsorshi 
of the National Bureau of Standards. The tables which have been made on 
able through the National Bureau of Standards are listed below. 

There is included in this list a publication on the hypergeometric and Legendre 
functions (MT15), prepared by the Bureau. 


MT1. Tasie or THE First Ten Powers OF THE INTEGERS From 1 To 1000: 


(1938) VIII-+-80 pages; heavy paper cover. 50 cents. 


MT2. Tastes OF THE EXPONENTIAL FUNCTION €?, 
The ranges and intervals of the argument and the number of decimal places in the entries 
are given below: 
Range of x Interval of x 

—2. 5000 to 1.0000 0. 0001 
1.0000 to 2.5000 - 0001 
2.500 to 5.000 - 001 
5.00 to 10.00 -01 


(1939) XV-+535 pages; bound in buckram, $2.00. 


MT3. Tastes or Circutar AND Hypersotic Sines AND Cosmes ror RADIAN ARGUMENTS: 


Contains 9 decimal place values of sin x, cos x, sinh x and cosh x for x (in radians) ranging from 
0 to 2 at intervals of 0.0001. 
(1939) XVII-+- 405 pages; bound in buckram, $2.00. 
MT4. Tastes or Smves AND Cosines FoR RADIAN ARGUMENTS: 


Contains 8 decimal place values of sines and cosines for radian arguments ranging from 0 to 25 
at intervals of 0.001. 
(1940) XXIX+-275 pages; bound in buckram, $2.00. 
MTS5, Tastes or Sine, Cosine, AND ExponentiAL IntTecRAts, Votume I: 
Values of these functions to 9 places of decimals from 0 to 2 at intervals of 0.0001. 
(1940) XXVI+- 444 pages; bound in buckram, $2.00. 
MT6. Taszes or Sinz, Cosine, AND ExponenTIAL InteGRALs, Votume II: 


Values of these functions to 9, 10, or 11 significant figures from 0 to 10 at intervals of 0.001 
with auxiliary tables. 
(1940) XXXVII-+-225 pages; bound in buckram, $2.00. 
MT7. Taste or Natura LocariruMs, VOLUME i: 
ithms of the integers from 1 to 50,000 to 16 places of decimals. 
(1941) XVIII+-501 pages; bound in buckram, $2.00. 
MT8. Tasres or Prosasiuitry Functions, Votume I: 


Values of these functions to 15 places of decimals from 0 to 1 at intervals of 0.0001 and from 
1 to 5.6 at intervals of 0.001. 
(1941) XXVIII+-302 pages; bound in buckram, $2.00, 


{Continued on p. 4 of cover} 
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MT9. Taste or Naturat Locaritums, Votume II: 


Logarithms of the integers from $0,000 to 100,000 to 16 places of decimals. 
(1941) XVIII+-501 pages; bound in buckram, $2.00. 


MTI10. Tasre or Natura Locarirums, Votume III: 


Logarithms of the decimal numbers from 0.0001 to 5.0000, to 16 places of decimals. 
(1941) XVIII+ 501 pages; bound in buckram, $2.00. 


MT11. Tasies or THE Moments oF Inertia AND Section Moput or Orpinary ANGLES, CHAN 
NEL8, AND Bus ANGLES wiTH Certain PLate ComsinaTions: 
(1941) XIII+ 197 pages; bound in green cloth. $2.00. 


MT12. Tasre or Natura Locarirums, Votume IV: 


Logarithms of the decimal numbers from 5.0000 to 10.0000, to 16 places of decimals. 
(1941) XXII+-506 pages; bound in buckram, $2.00. 


MT13. Taste or Sine AnD Costing INTEGRALS FOR ARGUMENTS From 10 To 100: 
(1942) XXXII-+-185 pages, bound in buckram, $2.00. 


MT14. Taszes or Prosasitry Functions, Votums II: 
Values of these functions to 15 places of decimals from 0 to 1 at intervals of 0.0001 and from 
1 to 7.8 at intervals of 0.001. 
(1942) XXI-+-344 pages; bound in buckram, $2.00. 
MT15. The h eometric and Legendre functions with applications to in’ uations of 
penne jan Chester Fo National Bureau of Seandards. Nadie 2 fous original 
ndwritten manuscript. 
(1942) VII+-319 pages; bound in heavy paper cover. $2.00. 
MT16. Tasre or Arc Tan X: 
Table of inverse tangents for positive values of the angle in radians. Second central differences 
are included for all entries. 


Range of x 


0 to 7 
7to 50 
50 to. 300 
300 to 2,000 
2,000 to 10,000 10 
(1942) XXV-+169 pages; bound in buckram, $2.00. 
MT17. Miscellaneous Physical Tables: 
Planck's radiation functions (Originally published in the Journal of the Optical Society of 
America, February 1940); and 
Electronic functions. 
(1941) VI+-58 pages, bound in buckram, $1.50. 


MT18. Table of the Zeros of the Legendre Polynomials of Order 1-16 and the Weight Coefficients 
for Gauss’s Mechanical Quadrature Formula: 
(Reprinted from Bull. Amer. Mathemical Society, October 1942.) 
5 pages. 25 cents. 

Payment is required in advance. Make remittance payable to the “National 
Bureau of Standards,” and send with order, using the blank form facing page 
3 of the cover. 

The prices are for delivery in the United States and its possessions and in 
countries extending the franking privilege. To other countries the price of 
MT1 is 65 cents; that of MT2 to MT16, inclusive, is $2.50 each; MT17, $1.75; 
MT18, 30 cents; remittance to be made payable in United States currency. 

Copies of these publications have been sent to various Government depositories 
throughout the country, such as public libraries in large cities, and colleges and 
universities, where they may be consulted. 

A mailing list is maintained for those who desire to receive announcements 
regarding new tables as they become available. 





